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Abstract. I discuss the spectrum of hadrons containing heavy quarks (b or c), and 
how well the experimental results are matched by theoretical ideas. Useful insights 
come from potential models and applications of Heavy Quark Symmetry and these can 
be compared with new numerical results from the ab initio methods of Lattice QCD. 

1 Introduction 

The fact that we cannot study free quarks but only their bound states makes 
the prediction of the hadron spectrum a key element in testing Quantum Chro- 
modynamics as a theory of the strong interactions. This test is by no means 
complete many years after QCD was first formulated. 

The 'everyday' hadrons making up the world around us contain only the 
light u and d quarks. In these lectures, however, I concentrate on the spectrum 
of hadrons containing the heavy quarks b and c (the top quark is too heavy to 
have a spectrum of bound states, see for example Quigg (1997a)) because in 
many ways this is better understood than the light hadron spectrum, both ex- 
perimentally and theoretically. The heavy hadrons only appear for a tiny fraction 
of a second in particle accelerators but they are just as important to our under- 
standing of fundamental interactions as light hadrons. In fact the phenomenology 
of heavy quark systems is becoming very useful; particularly that of B mesons. 
The study of B decays and mixing will lead in the next few years, we hope, to 
a complete determination of the elements of the Cabibbo-Kobayashi-Maskawa 
matrix to test our understanding of CP violation. CKM elements refer to weak 
decays from one quark flavour to another but the only measurable quantity is 
the decay rate for hadrons containing those quarks. To extract the CKM element 
from the experimental decay rate then requires theoretical predictions for the 
hadronic matrix element. We cannot expect to get these right if we have not pre- 
viously matched the somewhat simpler theoretical predictions for the spectrum 
to experiment. 

Here I will review the current situation for the spectrum of bound states 
with valence heavy quarks alone and bound states with valence heavy quarks 
and light quarks. The common thread is, of course, the presence of the heavy 
quark, but we will nevertheless find a very rich spectrum with plenty of variety 
in theoretical expectations and phenomenology. A lot of the recent theoretical 
progress has been made using the ab initio techniques of Lattice QCD. These 
are described elsewhere in this Volume ( Weingarten (1997)) along with recent 
results from Lattice QCD for the light hadron spectrum. 
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Quark model notation for the states in the meson spectrum will prove useful 
(baryons will not be discussed until section 3). The valence quark and anti-quark 
in the meson have total spin, S = or 1, and relative orbital angular momentum, 
L. The total angular momentum, which becomes the spin of the hadron, J = 
L + S. The meson state is then denoted by n 2S+1 Lj where n is the radial 
quantum number, n is conventionally given so that the first occurrence of that 
L is labelled by n—1 (i.e. n+1 is the number of radial nodes). L = is given the 
name S, L — 1, the name P, etc. To give J PC quantum numbers for the state 
(the only physical quantum numbers) we need the facts that P = (— 1) L+1 and, 
for C eigenstates, C = (— 1) L+S . In Table 1 a translation between 2S+1 Lj and 
J PC is provided. 



2S+l Lj 




'So 


o-+ 


3 Si 




'Pi 


1+- 


3 Po 


0++ 


3 Pi 


1++ 


3 P2 


2++ 



Table 1. J quantum numbers for quark model S and P states 



The ordering of levels that we see in the meson spectrum ( The Particle Data 
Group (1996)) is generally the naive one i.e. that for a given combination of quark 
and anti-quark adding orbital or spin momentum or radial excitation increases 
the mass. This is clearer for the heavy hadrons since, because of their masses 
and properties, the quark assignments arc unambiguous. For heavy hadrons it is 
also true, for reasons that I shall discuss, that the splittings between states of the 
same L but different S are smaller than the splittings between different values 
of L or n. To separate this fine structure from radial and orbital splittings it is 
convenient to distinguish spin splittings from spin-independent or spin-averaged 
splittings. Spin-averaged states are obtained by summing over masses of a given 
L and n, weighting by the total number of polarisations i.e (2 J +1). Examples 
are given below - they will be denoted by a bar. 

In Section 2 I begin with the phenomenology of mesons containing valence 
heavy quarks, the heavy-heavy spectrum. I shall discuss potential model ap- 
proaches to predicting this spectrum as well as more direct methods recently de- 
veloped in Lattice QCD. Section 3 will describe heavy-light mesons and baryons, 
both from the viewpoint of Heavy Quark Symmetry ideas and from Lattice QCD, 
using techniques successful in the heavy-heavy sector. Section 4 will give conclu- 
sions and the outlook for the future. 
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2 The Heavy-heavy Spectrum 

Figure 1 shows experimental results for bb and cc bound states (The Particle 
Data Group (1996)). They have been fitted on to the same plot by aligning 
the spin-average of the l 3 Po,i,2 states (x6 and Xc)- The spin-average \ mass is 
defined by 

M(x) = I [M( 3 P ) + 3M( 3 P0 + 5M( 3 P 2 )] (1) 

and this has been set to zero in both cases. The overall scale of bb meson masses 
is much larger than for cc but we see that this simply reflects the larger mass of 
the b quarks. The lightest vector state for bb is the Y produced in e + e~ collisions 
with a mass of 9.46 GeV. Its radial excitations are known as Y or Y(2S), Y or 
Y(3S) and so on. The radial excitations are separated from the ground state by 
several hundred MeV. For cc the lightest vector state is the J/ip or ip(lS) and 
this has a mass of 3.1 GeV. Its radial excitations are the tp' or tp(2S) and so on. 
Since the scale of Figure 1 spans 1 GeV it is clear that the splittings between 
states in both systems are very much smaller than the absolute masses of the 
mesons. 
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Fig. 1. The experimental heavy-heavy meson spectrum relative to the spin-average of 
the Xfc(l-P) and Xc(l-P) states (The Particle Data Group (1996)). 



It is also clear from Figure 1 that the radial excitations of the vector states 
in the two systems match each other very closely. In fact so closely that the 
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^(3770) which has vector quantum numbers cannot be fitted into a scheme of 
radial excitations of the ip system. It is thought to be not an S state but a D 
state (Rapidis et al (1977)). No bb D candidates have yet been seen). 

The matching of radial excitations is even better if we consider spin-averaged 
S states, 

M(S) = X - [MCSo) + 3M( 3 Si)] . (2) 

The 1 So state has only been seen for cc and is denoted r) c . The i] c (lS) lies below 
the vector state by 117 MeV and so the spin- average lies one quarter of this 
below the J/ip. As we shall see, this spin splitting in the bb system is expected 
to be much smaller. If we take a reasonable value for M(Y) — M{r)b) of 40-50 
MeV (see later), we would find the 15* levels on Figure 1 to be aligned to within 
10 MeV despite a difference in overall mass of a factor of 3. Similar arguments 
apply to the alignment of the 2S levels, although the agreement achieved there 
is not quite as good. 

The spin splittings within the Xb(lP) states (xbo, Xbi, Xb2) are much smaller 
than those within the Xc (\P) states, so that the spin splittings do depend on 
the heavy quark mass, vtlq, quite strongly. For example, we can take the ratio 
for IP levels: 

M( Xb 2) - M{ XW ) _ 53 MeV _ _ _. m 

M( Xc2 )-M( Xc0 ) - 141McV- U -' 8(i) - [6) 

Naively this looks very similar to the ratio of b and c quark masses if we take 
these to be approximately half the mass of the vector ground states, Y and 
J/ip. Then m c /m h w 0.33. This might imply a simple 1/roq dependence for spin 
splittings. However, the ratio between bb and cc does depend somewhat on the 
splitting being studied, indicating a more complicated picture. We have : 



M( X bi) - M( X bo) 
M(Xd) - M( Xc0 ) 

and 

M{ X b2) - M( X bi) 



0.34(2), (4) 



= 0.47(2). (5) 



M(xc2) - M(xci) 

The arrows shown on Figure 1 mark the minimum threshold for decay into 
heavy-light mesons, BB for bb and DD for cc. Three sets of S states and two sets 
of P states have been seen below this threshold for bb and two sets of S states 
and one set of P states for cc. Another set of P and two sets of D states are 
expected for bb (Eichtcn (1980), Kwong and Rosner (1988)). The states below 
threshold are very narrow since the Zweig-allowed decay to heavy-light states 
(see Fig. 2) is kinematically forbidden and they must decay by annihilation. This 
carries a penalty of powers of the strong coupling constant, ol s (Mq). These are 
then the states that we will concentrate on, because they can be treated as if 
they are stable (and none of the approaches which we will discuss allow them 
to decay). Vector states above threshold can still be seen in the e+e~ cross- 
section as bumps but they are much broader and the theoretical understanding 
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Fig. 2. Decay of a heavy-heavy meson to heavy-light mesons above threshold. 

of their masses requires a model for the inclusion of decay channels in the analysis 
(Eichten et al (1978), Ono et al (1986)). 

How can we understand the heavy-heavy spectrum below threshold? The 
fact that all the splittings are very much less than the masses, noted above, is 
critical. It implies that dynamical scales, such as the kinetic energy of the heavy 
quarks, are also very much less than the masses i.e. the quark velocities are 
non-relativistic, v 2 <C c 2 . Typical gluon momenta will be of the same order as 
typical quark momenta, rriQV. Thus typical gluon energies, toqw, are very much 
larger than typical quark kinetic energies, ttiqv 2 (Thacker and Lepage (1991)). 
The gluon interaction between heavy quarks will then appear 'instantaneous'. It 
can be modelled using a potential and energies found by solving Schrodinger's 
equation. In the extreme non-relativistic limit of very heavy quarks the spin 
splittings vanish. This was noticed above in the relation between bb and cc split- 
tings and will be discussed in more detail later. In this limit we need only a single 
spin-independent central potential to solve for the spin-averaged spectrum of 5 
and P states defined above. For a recent review of the history of the heavy quark 
potential see Quigg (1997b). 

2.1 The spin-independent heavy quark potential 

Perturbation theory for QCD gives a flavour-independent central potential based 
on 1-gluon exchange, which has a Coulomb-like form, 



where r is the radial separation between the two heavy quarks and a s is the 
strong coupling constant. The Coulomb potential cannot be the final answer 
because it would allow free quarks to escape. In addition it gives a spectrum 
incompatible with experiment in which the IP level is degenerate with 25. 

For a potential of the form V ~ r N with 2 > N > we have a IP level 
below 25, as we observe, and a ID level above 25 (as we see for charmonium). 



w 3 r 



(6) 



6 Christine Davics 



So, the addition of some positive power of r to the Coulomb potential can rescue 
the phenomenology (Grosse and Martin (1980)). The additional term is usually 
taken to be linear in r and thought of as a 'string-like' confining potential. This 
gives the simple Cornell potential of Eichten et al (1975): 

V(r) = -^+ar (7) 

with a called the 'string tension'. This can reproduce the observed spectrum 
reasonably well. Other successful forms for the heavy quark potential are the 
Richardson potential (Richardson (1979)): 

V(r) = J (8) 

in which a running strong coupling constant is included with non-perturbative 
behaviour at small q 2 (see also Buchmuller and Tye (1981)), and the Martin 
potential (Martin (1980), Grant, Rosner and Rynes (1993)): 

V(r) = Ar v with i/«0. (9) 

This last form, essentially a logarithmic potential (Quigg and Rosner (1977)), 
has no QCD motivation but is simply observed to work. All three potential 
forms can reproduce the bb and cc spin-averaged spectra reasonably well if the 
parameters arc chosen appropriately When this is done it is observed that the 
potentials themselves agree in the region r ~ 0.1 — 0.8 fm in which the \J (r 2 ) 
for the bound states sit (Buchmuller and Tye (1981)). 

It is interesting to compare the dependence of the energies of the states on 
the mass of the heavy quark, toq, in different potentials. This can be done 
for homogeneous polynomial-type potentials easily (see, for example Quigg and 
Rosner (1979), Close (1979)). Schrodinger's equation for the wavefunction \P is: 

h 2 V 2 1 

— — + V(r) \&{r) = E&(r). (10) 

E is the energy eigenvalue and /i, the reduced mass, toq/2 for the heavy- heavy 
case. For V = Ar N we can reproduce the same solution at different values of 
rriQ if we allow for a rescaling r — ► Ar. With this rescaling in place 

{-h 2 V 2 + A2fi\ N+2 r N } g'(Ar) = 2^X 2 E^(Xr). (11) 

The same solution (with rescaled r) will occur for different values of tuq if 

Aoc A1 - 1 /(2+^). (12) 

This gives a solution for E which varies as 

Eo,m Q N/{2+N) (13) 
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this tells us that 
and 



The values of E (and therefore splittings) will then be independent of tuq, as 
observed approximately, for N = 0. This corresponds to the Martin potential. 
The same result can be achieved by mixing N = 1 and N = -1 in the Cornell 
potential. Note that the Feynman-Hellmann Theorem guarantees that bound 
states fall deeper into the potential as the mass increases, dE/d^t < (Quigg 
and Rosner (1979)). For N > 0, E falls with toq; for N < 0, E increases in the 
negative direction. 

The Virial Theorem is helpful in extracting some dynamical parameters. It 
relates the mean kinetic energy to the expectation value of a derivative of the 
potential (see for example Quigg and Rosner (1979)): 

W = >f> d4) 

for homogeneous potentials. For N ~ i.e. V ~ logr we get K = a constant. 
Since 

2 

K 'k (15) 
(P 2 ) cx m (16) 

v 2 = 2{K)m2 (17) 

for a meson made of different quarks of masses mi and mi- v mi is the velocity 
of the quark of mass rm in the bound state. From fits using potential models a 
value of (K) of 0.37 GeV is found (Quigg and Rosner (1979)), giving 

v 2 

cin^p, — ~ 0.24 
c z 

v 2 

b in T, — - 0.07 

c l 

The quarks arc non-relativistic as we originally expected. For the logarithmic 
potential we also have the result that v 2 is independent of the radial excitation. 
For a Coulomb potential (p 2 ) decreases with increasing n, whereas for a linearly 
rising potential, (p 2 ) increases with increasing n. 

Potential model calculations of the bottomonium and charmonium spectra 
are reasonably successful. See Eichten and Quigg (1994) for a recent exam- 
ple, whose results are plotted in Figure 3. These include not only the central 
(Richardson) potential discussed here but also (perturbative) spin-dependent 
potentials to be described in section 2.2 to get spin splittings. In Eichten and 
Quigg (1994) parameters of the potential were fixed from a subset of states in 
the experimental cc spectrum. Typical deviations from experiment for the rest 
of the cc spectrum were 30 MeV; typical deviations in the bb spectrum were 
25 MeV. Since the b quark is significantly more non-relativistic in its bound 
states than the c quark one might expect to get better agreement for the bb 
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Fig. 3. The heavy-heavy meson spectrum from a recent Richardson potential model 
calculation (Eichten and Quigg (1994)). Circles and bursts show the calculated masses 
relative to the spin average of the Xb(lP) and Xc(l-P) states and the solid and dashed 
lines show experiment results, where they exist. 



spectrum using fitted parameters from that system. However, agreement for the 
cc spectrum would then be worse. In either case it is necessary to fit the param- 
eters of the phcnomcnological potential from some experimental information. 
Instead, the central potential V(r) can be extracted from first principles using 
the techniques of lattice QCD. 

In the iriQ — > oo limit the heavy quark is static. Its world line in space- 
time becomes a line of QCD gauge fields in the time direction. In Lattice QCD 
we break up space-time into a lattice of points and represent the gauge field 
by SU(3) matrices, U (Weingarten (1997), Montvay and Minister (1994)). The 
static quark propagator then becomes a string of U matrices, as in Figure 4. 



Fig. 4. The world line of a heavy quark on the lattice. 
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We can put a quark and antiquark together and join them up into a closed, 
and therefore gauge-invariant loop, called a Wilson loop (Figure 5). The value 
of the Wilson loop can be measured on sets of gauge fields {U} where a gauge 
field is defined on every link of the lattice. These are called configurations. The 
physically useful quantity is the matrix element of the Wilson loop between vac- 
uum states and this is obtained by averaging values of the Wilson loop over an 
ensemble of configurations where each configuration has been chosen as a typical 
snapshot of the vacuum of QCD. To obtain such an ensemble we must gener- 
ate configurations with a probability weighting e ~ SQCD and there are standard 
techniques to do this (Weingarten (1997), Montvay and Minister (1994)). 




Fig. 5. A Wilson loop. The Trace is over colour indices. 



The expectation value over such an ensemble of gauge fields of a Wilson loop 
of spatial size R is related to the heavy quark potential V(R). This is because 
the ensemble average is a Monte Carlo estimate of the path integral giving the 
matrix element of the operator which creates and destroys a static heavy quark 
pair at separation R on the lattice. The matrix element becomes exponentially 
related to the ground state energy of the quark anti-quark pair as the time extent 
of the Wilson loop, T, tends to infinity. Since R is fixed, and the quarks in this 
picture have no kinetic energy, this is simply the potential V(R) plus an additive 
self-energy contribution. 

„„., , f PC/Wilson Loop(C/)e~ s « CD 

(WfisonLoop) = JVUe -£j (18) 



[^ t (0) X t( J R)] t=o [v,(0) x ( J R)] t=T |0) 
|V(0)x(i?)|groi 
E = Vi att {R) = V(R) + constant. 



T ^°° |(0|V'(0)x(i?)|groundstate)| 2 e ET + higher order terms 
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How is the calculation done? Once the ensemble of gauge field configura- 
tions has been generated, Wilson loops of various different sizes in R and T are 
measured and average values of W(R, T) obtained. There is a statistical error 
associated with the number of configurations in the ensemble, i.e. how good 
an estimate of the path integral has been obtained. For a fixed R, W(R, T) is 
fitted to the exponential form above in the large T limit, extracting E. Away 
from T = oo higher order terms should be included in the fit which are ex- 
ponentials of excitations of the potential. There are a number of techniques to 
improve the values of E obtained, both the statistical error and any systematic 
error from fitting to an exponential form (see, for example, Bali, Schilling and 
Wachter (1997a)). Several of the techniques are similar to those used in direct 
calculations of the spectrum and are discussed in section 2.3. 

Once Viatt {R) is obtained it can either be used directly or a functional form in 
terms of R can be extracted to inform the continuum potential model approaches 
described above. The functional form usually used is that of the Cornell potential 
with e = 4a s /3 and an additive constant, V c : 

V latt (R) = aR-^ + V c . (19) 

The fit then yields the parameters a, e and V c . e is generally taken as a constant, 
although it is possible to determine the running coupling constant a s (R) from 
the short distance potential (UKQCD (1992b)). Often the running is mimicked 
by keeping e constant and adding an additional term, / /R 2 . This affects slightly 
the value of e obtained, as does the range of R included in the fit. A Martin form 
plus a constant, equation 9, does not fit Vi a tt (Bali, private communication). 

It is important to remember that Viatt, being obtained from the lattice, is 
measured in lattice units. To convert to dimensionful units of GeV we need to 
know the lattice spacing, a. This requires one piece of experimental information 
(see below). The separation, R, is also measured in lattice units, corresponding 
to a physical distance r — Rain fm. Thus the continuum potential V is obtained 
by 

V(r = Ra) = Vi at t(R)xa- 1 (20) 

However this expression should contain on the r.h.s. only the physical pieces of 
Viatt and not V c . 

V c is an unphysical constant which resets the zero of energy. It arises from 
corrections to the static quark self-energy induced by gluon loops sitting around 
the perimeter of the Wilson loop. These give a contribution to the logarithm of 
the Wilson loop proportional to its perimeter, V c (2R + 2T)/2. Thus the term V c 
appears as part of Viatt ■ In perturbation theory V c is a power series in the coupling 
constant a s , starting at 0(a s ), but is otherwise a constant in lattice units. From 
equation 20, its contribution to the continuum potential diverges on the approach 
to the continuum limit, a — > 0, and it should be subtracted from Viatt before 
equation 20 is applied. Another way to look at this is to notice that the heavy 
quark potential on the lattice is forced to zero at zero separation, Vj o tt(0) = 0, 
when the Wilson loop collapses to two lines on top of one another. Because the 
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U matrices are unitary we get (WilsonLoop) = 1 = e°. However, the continuum 
potential diverges in Coulomb fashion at zero separation so that V(0)a ^ 0. The 
physical pieces of the lattice potential will reproduce the continuum behaviour so 
to get Vi a tt(0) = will require an additive constant to shift the whole potential 
upwards. This is V c . 



> 
> 




Fig. 6. The heavy quark potential obtained from the lattice at two different values of 
the lattice spacing in the quenched approximation. The solid line is a fit of the form 19 
(Bali, Schilling and Wachter (1997a)). The potential and separation are given in units 
of the parameter ro (see text). 



Figure 6 shows recent results for the lattice potential plotted with the fit- 
ted form above, (19). The parameters extracted can be compared to those of 
phenomenological potentials. 

The Coulomb coefficient, e, is dimensionless and needs no multiplication by 
powers of the inverse lattice spacing to get a physical result, e is the coefficient 
of the 1 /R term but the discrete nature of the lattice changes 



1 4ttv-^ e l i- R 2 q. 



J = 2si„2i (21) 

<J#0 



Notice that this lattice form of 1/R is not rotationally invariant. At finite lattice 
spacing there are two alternatives. One is to fit this modified 'lattice' form of 
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1/R; the other is to correct for the discretisation errors in the naive lattice action, 
Sqcd which gave rise to them (Lepage (1996)). 

Most precision lattice calculations of the heavy quark potential have worked 
in the quenched approximation in which only the gluonic terms are included in 
Sqcd- This is equivalent to ignoring quark- antiquark pairs popping in and out 
of the vacuum (Weingarten (1997), Montvay and Minister (1994)). Recent calcu- 
lations (Bali, Schilling and Wachter (1997a), Bali and Schilling (1992), UKQCD 
(1992a)) have given a value of e around 0.3, which is rather smaller than the 
values that phenomenological potentials have used. For example, Eichten and 
Quigg (1994) use e = 0.54 in their Cornell potential fits. Part of this discrepancy 
can be traced to errors in the quenched approximation. When no qq pairs are 
available in the vacuum for screening, the strong coupling constant will run to 
zero too fast at small distances. Thus 

0t 8 {r)Q.A < &s(r) full theory 
V(r) Q . A . > V(r) full theory 

when V(r) is dominated by the Coulomb term. Calculations of the heavy quark 
potential that have been done on unquenched configurations (which usually con- 
tain two flavours of degenerate massive quarks in the vacuum, still not entirely 
simulating the real world), indicate that e is increased by about 10%. This gives 
a steeper potential at short distances, as in Figure 7. SESAM (1996) find eg. a. 
= 0.289(55) and e unquen ched = 0.321(100) without the use of the f/R 2 term 
in Equation 19. This doesn't then explain all of the difference between lattice 
values of e and phenomenological continuum values. 

Phenomenological potentials also implicitly include some relativistic correc- 
tions to the static picture that can be modelled simply as r-dependent additional 
potentials. The first such corrections contain a term inversely proportional to the 
square of the heavy quark mass multiplying a Coulomb potential, and therefore 
altering the effective value of e in an m-Q-dependent way. The coefficient of these 
corrections can be calculated on the lattice (Bali, Schilling and Wachter (1997a)). 
It is found that e becomes e + b/mq where b — (0.86(5)GcV) 2 , giving a signifi- 
cant increase (35%) to the effective value of e for charmonium but no change for 
bottomonium. This is illustrated in Figure 8, and supports the phenomenological 
use of different values for e in the two systems as a flavour-dependent dynamical 
effect. 

The string tension, a, describes the strength of the linearly rising part of the 
potential. It is dimensionful, so <Ji a tt = ca 2 . Using values of a from phenomeno- 
logical potentials ( Eichten and Quigg (1994) use ^fo w 0.43 GeV) allows us to 
fix a on the lattice and then convert all other dimensionful quantities to physical 
units. However, a and e are anti-correlated from the fitted form used in equation 
19 and this gives some bias. It is better to use instead the value ro obtained by 
setting r 2 F(r) to a fixed value. F(r) is the interquark force, obtained by differ- 
entiating the potential, and a suitable fixed value is 1.65 which corresponds to 
r « 0.5 fm (= 2.5GeV~ 1 when he = 1) (Sommer (1994)). Ensembles at different 
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Fig. 7. A comparison of the short-distance heavy quark potential obtained from the 
lattice on quenched and unquenched configurations by the SESAM collaboration. The 
potential and separation are given in units of the parameter ro (see text). Figure 
provided by Gunnar Bali. 



values of a are obtained by using different bare coupling constants in the action, 
Sqcd (Weingarten (1997), Montvay and Minister (1994)). However, the value of 
a for a given value of the bare coupling constant is not known a priori but has to 
be obtained by calculating a dimensionful parameter and comparing to experi- 
ment (or, in the case of a or ro above, to phenomenology). Fixing a is a critical 
step in a lattice calculation and introduces additional systematic and statistical 
errors into the quoted physical results. In the quenched approximation the value 
of a for a given gauge coupling, j3, will depend on the experimental quantity 
chosen to fix it, and so it is important to know what quantity was chosen when 
looking at lattice results. This point will be discussed further later. 

Given a value for a, Viatt — V c can be converted to GeV at separations, r, in 
fm (= GeV -1 ). This is then the physical heavy quark potential and it should be 
independent of the lattice spacing at which the calculation was done. Figure 6 
shows that this is true for current lattice results. 

Using the fitted form for the potential, the spectrum can be calculated by 
solving Schrodinger's equation in the continuum (Bali, Schilling and Wachter 
(1997a)). The heavy quark mass and the overall scale (given by the string ten- 
sion) need to be adjusted to optimise the fit to experiment. Including the rela- 
tivistic corrections to the potential described above and adjusting the value of 
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Fig. 8. A comparison of the heavy quark potential obtained from the lattice including 
the first relativistic corrections which yield an m^-dependent Coulomb term. The static 
(m,Q — * oo) potential is given by the solid line and those for b and c by dashed lines. 
Figure provided by Gunnar Bali, see Bali, Schilling and Wachter (1997a). 



e to mimic an unquenched result, yields average deviations from experiment of 
around 10 MeV for bottomonium and rather larger, as expected, 20 MeV for 
charmonium. The remaining systematic errors in the lattice potential (see sec- 
tion 2.2) could cause shifts of this size for bottomonium and make the 20 MeV 
deviations for charmonium look rather fortuitous. 



Exercise: Discuss how you would expect the potential appropriate to 
heavy baryons to behave. How would you calculate this on the lattice? 
(Thacker, Eichten and Sexton (1988)). 



2.2 The spin-dependent heavy quark potential 

As described above, the infinitely massive heavy quark is only a colour source; it 
carries no spin. To obtain spin splittings then we must move away from the static 
picture. A useful starting point is a non-relativistic expansion of the Dirac La- 
grangian which is appropriate for heavy quarks in heavy-heavy systems (Thacker 
and Lepage (1991)). This can be obtained by a Foldy-Wouthuysen-Tani trans- 
formation of the Dirac Lagrangian in Euclidean space (see for example Itzykson 
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and Zuber (1980)): 

C = tf(D t - -5- (22) 
2m Q 

CD 2 ) 2 ia 

- Ci ki +C2 84 (D ' E_E ' D) 

- C3 84' T - (DXE - EXD) - C4 2^' T - B --^ 

tp is a 2-component spinor with heavy quark and anti-quark decoupled. The 
mass term ip^niQip has been dropped. D is a covariant derivative coupling to the 
gluon field and E and B are chromo-electric and chromo-magnetic fields. The 
rest of the QCD Lagrangian for light quarks and gluons remains as usual. 

The terms in the Lagrangian can be ordered in powers of the squared velocity 
of the heavy quark using the following power counting rules for momentum and 
kinetic energy (Lepage et al (1992), Bodwin et al (1995)): 

D <~ p <~ uiqv 
K ~ mQV 2 



Then from the lowest order field equation 

] 

2toq 



(d t -igA i -^)^ = Q (23) 



we have 



gA 4 ~ dt ~ K = m,QV 2 
gE=[D t ,U]^ P K 



-ige ijk B k = [Di,Dj]~K 



2 2 4 

■ - m Q v . 



In C we then see that the leading order terms on the first line of equation 
22 are C(toqw 2 ) and these give spin-independent splittings in the heavyonium 
spectrum. On the second line are spin-independent terms of 0(m,QV 4 ) which are 
relativistic corrections to the leading terms. On the third line are spin-dependent 
terms also of (D(mQV 4 ). These are the leading terms as far as spin-splittings are 
concerned. So, as discussed earlier, spin splittings should be 0(v 2 ) times smaller 
than spin-independent splittings. This is equivalent to 1/wq behaviour, with a 
roughly constant kinetic energy, giving around 120 MeV for cc and 40 MeV for 
bb. Note that the er • (D x E) and a ■ B spin-dependent terms are of the same 
order because the chromo-magnetic field is suppressed by one power of v over 
the chromo-electric field in the power counting. 

Following these power-counting rules the number of operators to be included 
in C can be truncated at a fixed order in v 2 /c 2 and this is obviously a sen- 
sible thing to do if v 2 /c 2 <C 1. In describing heavy- heavy systems with this 
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Lagrangian, however, we have lost the renormalisibility of QCD. To obtain use- 
ful results we must put in a cut-off, A, to restrict momenta to p < A < vtiq. The 
excluded momenta e.g. in gluon loops will reappear as a renormalisation of the 
coefficients of the non-relativistic operators, the q of equation 22. The Cj can 
be calculated in perturbation theory (since they are dominated by ultra-violet 
scales for A >• Aq C d) by matching low energy scattering amplitudes of (22) to 
full QCD to some order in a s and p/rriQ. The Ci are all one at tree-level. 

I will describe two different, but related, approaches to the study of spin 
splittings in heavyonium. One is to take C of equation 22 and discretise it directly 
on the lattice - this is the NonRelativistic QCD approach (Thacker and Lepage 
(1991)). The second is to develop spin-dependent potentials from C to add into 
a Schrodinger equation, Hip = Eip, and solve for the splittings. 

ff = D|- c ^} + ^)+^ L ' S ^ S 2) ( 24 ) 

i=l 1 1 

where Vo(r) is the central potential from section 2.1 and V s a includes the spin- 
dependent potentials. Again one can take a phenomenological approach to the 
spin-dependent potentials, or extract them from the lattice. I will describe the 
potential approach first and then return to NRQCD in the next subsection. 





Fig. 9. The ratio of expectation values required for the spin-dependent potentials. F\ 
and F2 represent insertions of E or B as required for that potential. For some potentials 
these will be on the same side of the Wilson loop. The distance A should be large and 
t is summed over (see text). 



To extract spin-dependent potentials from QCD we start from the Wilson 
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loop which represents a static quark anti-quark pair at separation R (Eichtcn and 
Feinberg (1981), Pcskin (1983)). As discussed earlier, the heavy quark propagator 
in this case is simply a line in the time direction, from the simplest possible 
heavy quark Lagrangian, ip^D t ip. Imagine adding a perturbation er • B/2mg 
to the quark or anti-quark, such as would come from relativistic corrections to 
the propagator using the Lagrangian of equation 22. On one leg alone, zero is 
obtained by symmetry. If the perturbation is added to both legs and we sum 
over the time separations, t, between the two additions a new contribution to 
the potential is obtained of the form SiS 2 AV/mQimQ 2 . 



AV 



2 lim 

T^OO 



f 

Jo 



dt((B(O,0)B(R,t))) w 



(25) 



where (Q)w means the expectation value in the presence of the Wilson loop i.e. 
the ratio of the expectation value of the Wilson loop with the B field insertions 
to that without. This is easy to calculate using the methods of Lattice QCD 
(Michael and Rakow (1985), de Forcrand and Stack (1985)). Figure 9 illustrates 
this ratio for one value of t. An integration over t is required and this is approx- 
imated on the lattice by a sum (see Bali, Schilling and Wachter (1997a) for a 
recent description of the techniques used). The time separations of the insertion 
points from the ends of the Wilson loop, A, must be kept large to ensure that the 
spin-dependent contribution to the static propagation of a QQ pair is obtained 
in the ground state of the central potential; excited states must have time to 
decay away. 

The complete spin-dependent potential is given by (Eichtcn and Feinberg 
(1981), Chen, Kuang and Oakes (1995)): 



V sd = 



1 



Si 



2r \ m 2 - 



+ 



L 



Ql m Q2 

I ( Si + s 2 

r \m Q1 m Q2 
( Sr • rS 2 • r 



+ 



+ 



d V {r) + 2d 1 V 1 {r) 



(26) 



Ld 2 V 2 (r) 

1 Si-S 2 

\mQiniQ2r 2 3m Q1 mQ 2 



d 3 V 3 (r) 



Si • S 2 

3m Q1 m Q2 



d 4 V 4 (r) 




L do 



Ld 2 V 2 (r) 

y m Q1 mQ 2 ; 

The primes indicate differentiation with respect to the argument r of the 
different potentials. Note that the last two terms appear only for the unequal 
mass case. The di and di coefficients will be discussed below. In perturbation 
theory the di coefficients appear at 0(1), the di only at 0{a s ) and only for 
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m,Qi 7^ toq 2 . Vo is the central potential, discussed in section 2.1, and Vi, V~2, V3, 
V4 arc obtained on the lattice by calculating the following expectation values: 



^VI(R) = 2e ijk 



lim 



dtt 



E„ 



Bi 



/z- 



w 



(27) 



R 



V^iR) — £ijk hm / dtt 



/Z 



w 



(28) 



[RiRj ' Uij]V 3 (R) + UMR) = 2 hm f 



dt 



Bi 



/Z w (29) 



As in Figure 9, the denominator Zw is the expectation of the Wilson loop 
without E or B field insertions. 

A lattice discretisation of the E and B field strength operators is required. 
The simplest discretisation of F^ v (x) is to take the product of four U matrices 
around a 1 x 1 square in the /x, v plane starting from the corner x. This product 
is called a plaquette (Weingartcn (1997), Montvay and Minister (1994)); its 
hermitian conjugate should be subtracted and the resulting SU(3) matrix made 
traceless. Note that factors of g that would otherwise appear from equation 22 are 
absorbed into the lattice version of F M „. For the B field, a more symmetric version 
of this is to use, instead of one plaquette, the average of the four plaquettes 
around point x in the spatial plane perpendicular to B (see Figure 10). For 
E the spatial average of the two plaquettes at a given time is used. Sec Bali, 
Schilling and Wachtcr (1997a) for details. 

The central potential, as calculated on the lattice, is a spectral quantity, ap- 
pearing in the exponent of the exponential decay of a correlation function. It can 
therefore be directly interpreted as the continuum potential once converted to 
physical units. The spin-dependent potentials, in contrast, are calculated from 
the amplitudes of lattice correlation functions and undergo renormalisation when 
compared to continuum QCD. This renormalisation is visible in equation 26 as 
the di coefficients. These are functions of the Ci coefficients since the potentials 
are extracted by perturbing the Wilson loop with operators from equation 22 
(Chen, Kuang and Oakes (1995)). They reflect the matching required between 



The Heavy Hadron Spectrum 



1!) 



y 




Fig. 10. The sum of four untraced plaquettes around a point (the clover-leaf operator) 
that is used for a B z field insertion in a wilson loop for spin-dependent potentials (see 
text). 

this static/nonrelativistic effective theory and full QCD. In this case it is conve- 
nient to do the matching in two stages; full QCD to continuum effective theory 
and continuuum effective theory to lattice effective theory. 

For the first stage the a (and therefore di) have been calculated in leading 
order continuum perturbation theory (Eichten and Hill (1990), Falk, Grinstein 
and Luke (1991)). They depend logarithmically on the quark mass and the cut- 
off that is applied to the effective theory. The spin-dependent potentials also 
depend on the cut-off, but not m,Q, so that each term in V sc i becomes 

di(A,m Q )Vi(A). (30) 

We can use the di calculated in the continuum for the lattice calculation if we 
imagine the continuum effective theory at the same cut-off as the lattice cut-off 

(l/o)- 

For the second stage we then match between continuum and lattice effec- 
tive theories at the same cut-off. This provides an additional renormalisation 
which can be significant because of the non-linear relationship between the con- 
tinuum and lattice gauge holds (Lepage and Mackenzie (1993)). This gives rise 
to additional tadpole diagrams in lattice perturbation theory. They have a uni- 
versal nature and can be thought of (even beyond perturbation theory) as a 
constant factor multiplying each gauge link. Equivalcntly the renormalisation 
can be viewed as arising from the additional perimeter self-energy contributions 
when the E and B field insertion arc in place (de Forcrand and Stack (1985)). 
A method to take account of the renormalisation directly on the lattice involves 
calculating, instead of the ratio in Figure 9, the product of ratios in Figure 11 
(Huntley and Michael (1987)). The additional perimeter/tadpole terms from the 
insertions are thereby cancelled out, and it is hoped that any remaining lattice 
renormalisation is negligible. 

Once the spin-dependent potentials are calculated from Figure 11 and mul- 
tiplied by the appropriate di, they can be inserted into a Schrodinger equation 
and the energy shifts from the spin-independent states can be calculated (Bali, 
Schilling and Wachter (1997a)). They depend on the functional form of the spin- 
dependent potentials and on the expectation value of the spin and orbital angular 
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Fig. 11. The ratio of expectation values used for the spin-dependent potentials, taking 
account of renormalisation required to match to the continuum (Huntley and Michael 
(1987)). F\ and F 2 represent insertions of E or B as required for that potential. 



momentum operators of equation 26 for a particular state. To calculate these 
the following equations are useful (for the last relation see Kwong and Rosner 
(1988)): 



(Si • S 2 ) 


1 

~ 2 


S{S + l)-\ 






(L-SO 


= <L 


■S 2 ) = 1(L-S) 




(L-S) 


= \[J{J+l)-L{L + l) 


-5(5 + 1)] 


(Sij) 


= 4<3(Sj ■ nKS,- ■ n) - S< 






= 2<3(S ■ n)(S ■ n) - 


-S 2 > 





[12(L • S) 2 + 6(L • S) - 45(5 + 1)L(L + 1)] 
~ (2L- l)(2L + 3) 

The results for the expectation values are tabulated for 5 and P states in 
Table 2. It is clear from this table that the only potential contributing to the 
hyperfme splitting between the 3 5i and 1 S states (M(T) - M(rj b ), M{J/&) - 
M(rj c )) is Vi. For the splittings between P states, all the spin-spin and spin-orbit 
potentials can contribute in principle. Notice how the spin-averaging described 
at the beginning of section 2 removes all the spin-dependent pieces, to obtain the 
spin-independent spectrum. To remove V4 terms from P states the spin-average 
must be taken including the l P\. 
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'So 


3 Si 


'Pi 


3 Po 


3 Pi 


3 P 2 


(Si ■ S 2 ) 


3 
4 


1 

4 


3 
4 


1 

4 


i 

4 


i 

4 


(L-S) 











-2 


-1 


1 


{Sa) 











-4 


2 


2 
5 



Table 2. Expectation values for combinations of spin and orbital angular momentum 
operators needed for spin splittings in heavy-heavy bound states. 



What functional form do wc expect for the different spin-dependent poten- 
tials? In leading order perturbation theory (one gluon exchange): 

V = -C F — (32) 
r 

V x = 

V 2 = C ^ 
V * = 3C ^ 

V 4 = 8irC F a s 6( 3 \r), 

with Cp = 4/3. The 'same-side' (see equation 27) spin-orbit interaction, V\, is 
absent; the 'opposite-side', V 2 is simply V . The form of V4 implies that it is only 
effective for states with a wavefunction at the origin i.e. S states. It gives for the 
3 Si - 1 S splitting, 

l^fMO)!'. (33, 

We do not then expect any splitting induced by the V4 term between the 1 P\ 
and 3 Pi states, so the 1 P\ mass should be at the spin-average of the 3 P states. 
The following inter-relationships between potentials are also useful. 



V 2 -V 1 = V (34) 

V 3 (r) = ^--V 2 '(r) (35) 
V 4 (r) = 2W 2 V 2 (r) (36) 



Equation 34 is the Gromes relation (Gromes (1984)), derived from Lorentz in- 
variance and as such always true. Equations 35 and 36 hold for any vector-like 
exchange (such as single gluon) but only to leading order; they do not survive 
renormalisation of the potentials when the cut-off on the effective Lagrangian 
of equation 22 is changed. In particular V\ and V 2 then mix (Chen, Kuang and 
Oakes (1995)). 

A crucial ingredient missed in the perturbative analysis is the confining part 
of the central potential, Vq, and this can reappear in the Vi. The nature of this 
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confining term is important. General considerations (Gromes (1977)) show that 
it can only arise from vector and/or scalar exchange, but these two possibilities 
yield quite different accompanying spin-dependent potentials. A vector exchange 
gives rise to V 2} V 3 and V4, a scalar exchange only to V\ (Gromes (1988)). In 
both cases, a constant term a appears in V 2 ' — V{ from the Gromes relation. 
A useful quantity to study in this respect is the ratio of p state splittings: 

M( 3 P 2 ) - M( 3 PQ 
' M( 3 P!)-M( 3 Po)' 

Experimentally this ratio takes the value 0.48(1) for cc(\P) and 0.66(2) for 
bb(\P) and 0.58(3) for bb(2P). For pure L- S interactions p is simply related to a 
combination of expectation values of L • S since, considering the spin-dependent 
potentials as a perturbation on the spin-independent one, the expectation value 
of Vi in all the P states is the same. This then gives p = 2. Similarly a pure tensor 
V3 interaction gives p — -0.4. These are clearly wrong; we require a mixture of 
spin-orbit and tensor terms. For the leading order perturbative potentials in 
equation 32 we can also calculate p exactly because all the expectation values of 
Vi reduce to cancelling terms of the form (r~ 3 ). This gives p = 0.8, larger than 
all the experimental values. The confining term should then appear in such a 
way as to reduce p. 

This is possible if we make the assumption that the confining term grows 
linearly with r as or; such a rapid rise implies a scalar exchange (Gromes (1988)). 
V\ = —or and V2 = —CF&s/r will satisfy the Gromes relation, p becomes 

= 1 8a s (r- 3 ) -5/2a(r- 1 ) 
9 5 2a s (r- 3 ) - l/4a(r" 1 ) 1 ' 

and for positive expectation values this will be less than 0.8, in agreement with 
experiment (Henriqucs, Kcllctt and Moorhouse (1976)). The a term will be more 
effective for longer range wavefunctions such as cc and bb(2P) giving a smaller 
value of p than for bb(lP). A vector confining potential would lead to the term 
proportional to a appearing in Vi with opposite sign as well as additional V3 
terms, so that p > 0.8 (Schnitzcr (1975)). Of course this docs not rule out a 
mixture of long-range vector and scalar pieces. 

The lattice calculation of the spin-dependent potentials confirm the behaviour 
above explicitly, and show (within errors) that the long range confining potential 
is purely scalar (Huntley and Michael (1987)). V3 and V4 are found to be very 
short-range with V3 showing 1 / R? behaviour and V4 approximating a 8 function 
on the lattice. V{ is approximately constant at the value —a taken from the 
central potential, whereas V 2 — > at large R. In addition V x has a small attrac- 
tive 1/R 2 piece (see Figure 12 from Bali, Schilling and Wachter (1997a)). which 
arises from the mixing between V\ and V 2 and its size changes as the lattice 
cut-off (1/a) changes, along with the Coulombic 1/R 2 term present in V 2 . 

There is no exact Gromes relation on the lattice (Bali, Schilling and Wachter 
(1997b)), but it should be restored in the continuum limit. This relation does 
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Fig. 12. The spin-orbit potential, — Vi , at two different values of the lattice spacing 
together with a fit curve of the form a + h/R 2 . (Bali, Schilling and Wachter (1997a)). 



in fact work well on the lattice at current values of the lattice spacing and 
this is a non-trivial check of the lattice renormalisation procedure of Figure 1 1 
(Huntley and Michael (1987)). This renormalisation is done for the left hand 
side of equation 34 but not for the right. See Figure 13 from Bali, Schilling and 
Wachter (1997a). 

As discussed earlier, the spectrum from this lattice potential yields deviations 
at the 10 MeV level for bottomonium and the 20 MeV level for charmonium. 
Systematic errors in the charmonium case are rather larger than this, however. 
The di coefficients have large perturbative corrections for the lattice cut-off used 
(m c a < 1) and so large uncertainties. Ci is set to 1 in equation 24; unknown 
perturbative corrections to that coefficient could induce 50MeV shifts in the 
charmonium spectrum (Bali, Schilling and Wachter (1997a)). As mentioned in 
section 2.1 there are also relativistic corrections to the spin- independent central 
potential (Barchielli, Brambilla and Prosperi (1990)). These can be calculated 
from expectation values of Wilson loops with E and B insertions in a similar way 
to that for the spin-dependent potentials above. The results modify the central 
potential for charmonium quite strongly, again indicating that unknown higher 
order corrections could be significant for that system. 

To go beyond the corrections discussed here in the potential model approach 
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Fig. 13. The difference of spin-orbit potentials, V 2 and V 1 on the lattice, compared 
to the expectation from the central potential according to the Gromes relation. (Bali, 
Schilling and Wachter (1997a)). 

is hard; higher order insertions into Wilson loops cannot be reduced to the 
form of an instantaneous potential. We need instead more direct methods of 
calculating the spectrum. This must be done on the lattice and will be discussed 
in the next subsection. 

Exercise: Fill out Table 2 to include D states. 

2.3 Direct measurement of the bottomonium spectrum on the 
lattice 

A direct calculation of the heavyonium spectrum on the lattice at first sight 
seems rather hard. There is a large range of scales in the problem, all the way 
from the heavy quark mass to kinetic energies within bound states (w Aqcd)- 
To cover these properly in a lattice simulation would require a -1 ^> tuq and the 
number of lattice points on a side, L » itiq/Aqcd- 

As we have seen from the previous sections, however, the quark mass itself 
is not a dynamical scale, simply an overall energy shift. We only actually need 
to simulate accurately the important scales for the bound state splittings, pq 
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and K. This leads us to work with a lattice with a -1 < O(mg) and make use 
of the non-relativistic effective theory of equation 22. This Lagrangian can be 
discretised on the lattice (Thacker and Lepage (1991), Lepage et al (1992)) and 
applied using similar techniques to those for handling light quarks on the lattice 
(Weingarten (1997), Montvay and Minister (1994)). Details will be discussed 
below. 

There is an important difference between the NRQCD approach and the 
potential model approach of the section 2.2. That approach starts from the 
static theory and so can only produce the potential; the missing kinetic energy 
terms are of equal weight (in powers of v 2 /c 2 ) in the spectrum and they must 
be added in subsequently in a Schrodinger equation. The NRQCD calculations, 
even at lowest order, include both the tp^ D t ip term and the tp^ D 2 /2niQip terms 
and yield the spectrum directly; the existence of a potential is not invoked at any 
stage. This means that the NRQCD approach can be fully matched to QCD and 
handle the sub-leading effects from soft-gluon radiation that eventually cause a 
potential model picture to break down through infra-red (long time) divergences 
(Appelquist et al (1978), Thacker and Lepage (1991)). We will find potential 
models useful for guiding NRQCD calculations, nevertheless. 

The NRQCD approach uses the Lagrangian of equation 22 as an effective 
theory on the lattice (Lepage et al (1992)). It can reproduce the low energy 
(p < 1 1 a) behaviour of QCD, but the couplings, Cj, must be adjusted from their 
tree level values of 1 to compensate for neglected high momentum interactions. 
In principle this can be done in perturbation theory by matching scattering 
amplitudes between lattice NRQCD and full QCD in the continuum (here a one- 
stage matching is used). The Ci will have an expansion in terms of a s (l/a). They 
will differ from the Ci of the static approach discussed earlier since the p 2 /m,Q 
term in the heavy quark propagator will give additional explicit l/m,Qa terms 
which diverge as a — ► 0. In this way it is clear that we cannot take a continuum 
limit in NRQCD; we can only demonstrate that results are independent of the 
lattice spacing at non-zero lattice spacing. This is sufficient for them to make 
physical sense, and to be compared to experiment. 

One problem for lattice NRQCD is the possible large renormalisations, Cj, 
which come from tadpole diagrams. This was discussed earlier in connection 
with the renormalisation of the spin-dependent potentials in the static case. The 
tadpoles appear with every occurrence of a gluon link field and can be taken 
care of by renormalising each gauge link by a factor u as it is read in, 



and then using the the renormalised gauge link everywhere instead of the origi- 
nal. 

uq represents how far the gluon links are from their continuum expectation 
value of 1. The easiest quantity to use to set uo is the plaquctte. Since it contains 
four links we have: , 




(39) 




(40) 
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A possibly better motivated value is that in which we look at a single link field 
and maximise its value by gauge-fixing. This should be most effective at isolat- 
ing (by minimising) the true gauge-independent tadpole contribution (Lepage 
(1997)). The gauge in which this happens is Landau gauge: 

"0 = U 0L = -Tr(f7 M ) Landau gauge- (41) 

The difference between the two w s is small in lattice perturbation theory. Mea- 
sured (non-perturbative) values on the lattice differ by a few percent at moderate 
values of the lattice spacing. 

Once we have renormalised the gauge fields to take account of the tadpole 
contributions ('tadpole-improvement') we would hope that remaining corrections 
to the Ci are small. They have been calculated in perturbation theory for those 
terms which contribute to the heavy quark self-energy (Morningstar (1994)). 
The dispersion relation for the heavy quark is required to be 

with A an energy shift and m r the renormalised quark mass, and this fixes the 
coefficient c\ in the lattice discretised version of equation 22. Figure 14 shows 
that the 0(a s ) coefficient of c\ is small, its magnitude less than 1, until mga 
is less than about 0.8, when it starts to diverge. This is a sign of the power 
ultra-violet divergences of NRQCD mentioned above; we must stay at values of 
a where mga > 0.8. Without tadpole-improvement the 0(a s ) coefficients are 
all much larger than 1 (Morningstar (1994)), showing that tadpole-improvement 
has captured most of the renormalisation. The results I shall describe here use 
tadpole- improvement and all Ci then set to 1. 

The NRQCD Lagrangian is discretised on the lattice in the standard way 
(Wcingartcn (1997), Montvay and Minister (1994)). Derivatives are replaced 
by finite differences, and E and B fields by clover terms. In the process, all 
appearances of tuq are replaced by the bare quark mass in lattice units, niQa, 
and powers of g are absorbed into the lattice fields. The lowest order terms in 
the Lagrangian density (in lattice units) become 

D t i/H -> Ut^t+i ~ A (43) 

d 2 ; E^A+i + t^-r^ 

2toq 2m,Qa 

Each Ufj. field here is understood to have been divided by uo already. 

Then the calculation of the heavy quark propagator is very simple. The 
propagator as a function of spatial indices on a given time slice is related to that 
on previous time slice by an evolution equation : 



U t G t+1 -G t = -aHG t (44) 



G t+1 = U}(l-aH)G t 
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Fig. 14. The C(a s ) coefficients of various terms in the NRQCD Lagrangian, calcuiated 
in lattice perturbation theory after tadpole improvement. A corresponds to the energy 



shift and to the D 4 term. B corresponds to the mass renormalisation and c£' to 
the Df term, here denoted c 5 . The vertical lines represent discontinuities when the 
value of the stability parameter, n, is changed. (Morningstar (1994)). 



where aH is the Hamiltonian, for example, the lowest order D 2 /2mQ term, dis- 
cretised on the lattice as in equation 43. This enables the heavy quark propagator 
to be calculated on one pass through the lattice in the time direction starting 
with some source for the propagator on time slice 1. This simplicity can be traced 
back to the simple first order time derivative in equation 22; calculations of rel- 
ativistic quark propagators take many sweeps through the lattice (Weingarten 
(1997), Montvay and Minister (1994)). 

One technical problem with (44) is that it can become unstable for modes 
for which H approaches 1. In the free case for the lowest order Hamiltonian we 
have 

3E J 4sin 2 (p l a/2) 



H 



(45) 



2niQa 

where i runs over the Fourier modes and the maximum value for momenta close 
to the lattice cut-off is 6/rriQa. This would limit values of rriQa to be greater than 
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6. Instead we can stabilise the evolution by adding terms which have an effect at 
the cut-off scale but are not important for the physically relevant momenta well 
below the cut-off. This gives rise to an evolution equation (Thacker and Lepage 
(1991)): 

and the stability requirement is now rriQa > 3/n so reasonable values of rag a 
of 0(1) can be reached for suitable n. In fact it is only important to stabilise 
the lowest order term Ho like this; the higher order terms of the Lagrangian 
of equation 22 can be added in straightforwardly. For example, we can write 
(Lepage et al (1992)): 



2n J 1 V 2n , 
with aSH the lattice discretisation of : 

SH = -ci^) +c 2 ^ r (D-E-E-D)-c 3 -^ r cr-(DxE-ExD)-C4-^cr.B. 

8toq 8toq 8mg 2toq 

(48) 

Another technical problem which faces all lattice calculations is that of dis- 
cretisation errors. These arise from the use of finite differences on the lattice 
to approximate continuum derivatives. They mean that even physical lattice 
results, expressed in GeV for example, will depend upon the lattice spacing. 
The dependence will be as some power of a typical momentum scale in lattice 
units. Since the momenta inside heavyonium systems are quite large, 0(1 GeV), 
discretisation errors will cause a problem if they are not corrected for. 

This is achieved by improving the discretisation of derivatives to include 
higher order terms. For example, ignoring gauge fields, 

a 2 Dl latt iP x = <P x+i + - 2^ (49) 
= (e aD i,com - 1) (1 - e -«A,co„t) ^ 

[---] 
L 6 4 J 



= ( o 2 D? t + a 4 Dl cont £ - 1] . . . ) V*. 



giving 0(a 2 ) errors relative to the leading term. A better discretisation is then 

a 2 Dl latt = a 2 Dl latt - ±a*Bf tlatt (50) 

where a 2 D 2 latt is given by the naive finite difference. a 2 D 2 latt has errors at 
relative 0(a 4 ). 

The other operator that appears in the leading order terms is the time deriva- 
tive operator, D t . Any correction to D t that looks like D\ would upset our simple 
evolution equation. Instead the way to correct D t is to require that the time-step 
operator be 

G t +i - e- aH G t . (51) 
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In fact the modified evolution equation (46) is closer to this than (44), and 
would be correct for the kinetic terms in the n — > oo limit. The gauge potential 
will appear automatically exponentiated from the appearance of l/J . The only 
correction that then needs to be made at the next order is to correct for a Hq jn 
term that appears when (46) is expanded out and compared to (51) for H = Hq. 
This correction can be made by replacing H by 

H a = H -^-. (52) 
An 

The discretisation corrections discussed here, when added in to 5H look a lot 
like relativistic corrections. We can apply the same power counting arguments 
as before to get an idea of their relative size. From (52) we will have a correction 
of 0(amQV 4 ). Relative to H this is 0(arriQV 2 ), and for arriQ <~ 1 this is 0(v 2 ), 
the same as the relativistic corrections. Similarly for the term from (50). Thus 
it is only sensible to correct for the discretisation corrections up to an order 
comparable with the order of relativistic corrections being included. For the La- 
grangian of 22 which has the first spin-independent relativistic corrections we 
need only the first spin-independent discretisation corrections described above. 
It might be true on coarse lattices, with toqo > 1, that higher order discretisa- 
tion corrections should be kept. This can be decided by using potential model 
expectation values to better estimate their size (Lepage (1992)). 

There are additional 0(a 2 ) errors from the gluon fields that appear in all the 
covariant derivatives coupling to the heavy quarks, if the gluon fields have been 
generated using the standard Wilson gluon action (Weingarten (1997), Montvay 
and Minister (1994)). These errors can be treated perturbatively and corrected 
for at the end of the calculation provided they are small (Davies et al (1995a)). 

The coefficients of the additional terms introduced by the discretisation cor- 
rections in (50) and (52) must again be matched to full QCD. As before, they 
should be tadpole-improved to remove the largest part of the renormalisation 
of lattice NRQCD, and remaining renormalisations can be calculated in lattice 
perturbation theory. The improvement from (52) can be added directly to the 
existing relativistic correction to give the operator 



"fc^^f fl + ^)- (53) 



The 0(a s ) corrections to c\ were discussed above and are shown in Figure 14. 
The improvement from (50) gives 

5H,=cM — . (54) 

2Am Q a 

The 0(a s ) corrections to C5 are also shown in Figure 14 and they are confirmed 
to be small after tadpole- improvement (see Morningstar (1994) and note that c 5 
is there called c 2 ). 



30 Christine Davies 



Once the lattice NRQCD Lagrangian (including discretisation corrections) 
has been chosen and the quark propagator Gt calculated for a given gluon field 
configuration, then Gt and the anti-quark propagator G t can be put together to 
make mesons. This procedure is identical to that used in lattice calculations of 
the light hadron spectrum. The only difference is that the meson operator : 

^(x 1 )^(x 1 -x 2 ) X tA (x 2 ) (55) 

has a spin part, fi, which is only a 2 x 2 matrix, rather than the relativistic 
4 x 4. Q is the unit matrix for 5 = mesons and a Pauli matrix for S = 1. In 
addition we have a much better idea from potential models of what form the 
spatial operator should take, than we do in light hadron calculations. This will 
be discussed below, ip> and x^ are the quark and anti-quark creation operators 
respectively, matched in colour, denoted A, for a colour singlet. 

Then the meson correlation function is calculated as an average over the 
ensemble of gluon field configurations (Weingarten (1997), Montvay and Miinster 
(1994)): 

((x^VOtW^X^o) = {Tr[GQ^^G^Q]) (56) 
T => 00 <P ie - BlT + <P 2e -^ T + .... 

Ei and E 2 are the energies of states in lattice units, E p hy S a. E\ is the energy of 
the ground state in that Q, <p channel, and E 2 is the energy of the first radial 
excitation etc. We can project onto different meson momenta at the annihilation 
time point, T, to obtain the dispersion relation, £ as a function of p (Davies et 
al (1994a)). 

Because it is very important in heavyonium physics to calculate radial excita- 
tion energies, we need to optimise the calculation of E\ and E 2 . The coefficients 
<Pi and <P 2 in (56) represent the overlap of the mesonic operator used with that 
state, <Pi = (0\^ Q(f>x^ \i) , see equation 18. We can then adjust <Pi and <P 2 by 
changing (f> in the mesonic operator. For S states we could use an operator in 
which both tp and x appear at the same point but this would have overlap with 
all possible excitations and a poor convergence to the ground state. Instead we 
must tp and x at separated points with <fr a 'smearing function'. Potential model 
wavefunctions represent a good first approximation to the spatial distribution of 
quark and anti-quark in heavyonium (Davies et al (1994a)). To make use of these 
wavefunctions on the lattice (specifically to set </> equal to the wavefunction) re- 
quires us to fix a gauge otherwise the meson operator will not be gauge-invariant 
and will vanish in the ensemble average. The best gauge to use is Coulomb gauge 
since this (being the 3-dimensional version of the lattice Landau gauge discussed 
earlier) is the gauge in which the spatial gluon field is minimised, and the co- 
variant squared spatial derivative most like the Schrodingcr p 2 . We then gauge 
transform the lattice gluon fields to Coulomb gauge and use different <p as sim- 
ple functions of spatial separation for different radial and orbital excitations. For 
the ground state </> should maximise <Pi and minimise <P 2 , <?3 etc, and for excited 
states <p should maximise $ 2 and minimise <Z>3 etc. For each <j> a new quark 
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propagator must be calculated in this approach since the fastest way to make 
the meson operator at the initial time slice is to use (j)(x.) as a source for the 
evolution equation for G t - This G t is then combined with a G\ in which a delta 
function at the spatial origin was used as the source. The same source </> can be 
used for the different spin states (to the extent that spin-dependent effects on 
the wavefunction are relativistic corrections and therefore small) and the factors 
of ft at the initial time slice inserted as Gt and G\ are being combined. In this 
way all 2S+1 Lj states can be made with as many radial excitations as required 
(Davies et al (1994a)). 

It is important to realise that the results for E\ and Ei are not affected by 
the choice of <f>] they can simply be obtained more efficiently by good choices. 
Methods other than that above have also been used; these include building meson 
operators out of a quark and anti-quark joined by a string of gauge fields in a 
gauge-invariant way (Manke et al (1997)); and calculating the propagators from 
delta function sources at a number of spatial points at the initial time and 
working out the optimal cf) at the end using a variational method (Draper et al 
(1995)). 

However good the choice of <f> 1 each meson correlation function will contain 
several exponentials and a multi-exponential fit must be performed to extract 
them. This is described with technical details in Davies et al (1994a). In gen- 
eral the nth exponential is obtained reliably from an n + 1-exponential fit. In a 
potential model approach to the spectrum, using orthogonal wavefunctions, it 
is easy to get very precise results for radially excited states. In lattice NRQCD 
it is much harder because the ground state will take over exponentially if it is 
present at all in an excited meson correlation function. In addition the variance 
of such a correlation function will be dominated by the ground state so that the 
ratio of the signal for the excited state compared to noise will fall exponentially 
(see, for example Lepage (1989)). 

The fits to the zero momentum meson correlation functions yield a very accu- 
rate set of energies in lattice units but these cannot be immediately converted to 
absolute energies (although splittings can) because the zero of energy has been 
reset by the absence of the mass term in equation 22. To calculate the spectrum 
we must shift all the lattice energies by a constant and then convert to physical 
units by multiplication with a -1 . 

Determining the lattice spacing is actually easier in heavyonium than for 
light hadrons. We can make use of the fact, stated before, that the radial and 
orbital splittings are independent to a very good approximation of the heavy 
quark mass. This means that we can use one of these splittings, e.g. the IP — IS 
splitting, to determine a -1 , without having necessarily tuned our heavy quark 
mass very well. In the absence of an experimentally determined spin-average S 
state mass for bottomonium we set 

a' 1 = — ^rGeV. (57) 

aE(xE) - aE{T) V ' 

where the denominator is the difference between the lattice energies at zero 
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momentum of the spin average of the ground Xb states and the T. Given a 1 we 
can now convert all differences in energy to splittings to GeV. We can also use 
E{T) -E{Y) to set a" 1 . 

To tune the bare lattice heavy quark mass, mg, to the appropriate value for 
the b quark we study the dispersion relation for mesons at finite and small mo- 
menta, where the heavy mesons are non-relativistic. The absolute meson mass 
(e.g. for the T) is given not by the energy at zero momentum but by the denom- 
inator of the kinetic energy term : 

2 2 

aE r {p) = aE r (0) + + . . . . (58) 

Higher order relativistic corrections can also be added to this formula. We adjust 
vfiQd in the Lagrangian until the T mass comes out at 9.46 GeV within statistical 
errors, using the a -1 determined from the splitting above. Now it is clear that if 
the splittings used for determining a -1 did depend strongly on m,Qa this would 
be a tricky iterative procedure. It would require complete calculations at several 
different values of rriQa, as is generally undertaken in light hadron calculations. 

This procedure gives us also the shift of the zero of energy, aMr — aEr(0). 
It should be independent of the meson studied and so, once calculated, can be 
applied to all mesons. That is, when divided by 2, it can be applied as a shift per 
quark. This is what allows us to convert differences in zero momentum energies 
on the lattice directly to splittings in physical units, given a -1 . 

The shift obtained can be compared to that calculated in lattice perturbation 
theory from the heavy quark self-energy. The energy shift is given in lattice units 
by 

2{Z m m Q a - E a) (59) 

where Z m is the mass renormalisation. Z m and E^a are given by perturbative 
expansions, Z m = 1 + a s B + . . . and E a — a s A + . . .. Again it is clear that 
A and B are smaller when a tadpole-improved lattice Lagrangian is used, (see 
Morningstar (1994) and Figure 14). The shifts obtained on the lattice agree well 
with the perturbative estimates (Davies et al (1994a), Davies (1997)) when a 
physical scheme for the lattice coupling constant is used (Lepage and Mackenzie 
(1993)) and allowance is made for unknown higher order terms. See Table 3. 

Note that if we take m,Qa to oo in this calculation aE will become V c /2 where 
V c is the unphysical self-energy part of the heavy quark potential discussed ear- 
lier. Again agreement between perturbation theory (Morningstar (1994), Duncan 
et al (1995)) and potential model results (Bali and Schilling (1992)) is reasonable 
given that aE starts at 0(a 8 ) and is only known to this order. The effect of 
tadpole-improvement is easy to work out in this case. If the Wilson loops were 
calculated with tadpole-improvement (not usually done) then each would be 
divided by uq and the loop would pick up a factor («o)~ 2T from links in the time 
direction. Then V c — > V c + 2\nu . Thus to compare the perturbative value of 
aE (m,Qa — ► oo) calculated with tadpole-improvement to the non-perturbative 
values of V c /2 calculated without tadpole-improvement we must subtract lnu 
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p 


rriQa 


Perturbative shift 


Non-perturbative shift 


5.7 


3.15 


7.0(6) 


6.54(7) 


6.0 


1.71 


3.5(2) 


3.49(3) 


6.2 


1.22 


2.5(2) 


2.58(3) 



Table 3. Energy shifts for a heavy quark in lattice NRQCD. Results are given for the 
non-perturbative lattice calculation of aMr — aEr(0) and for the perturbative shift 
of equation 59 for three different values of the lattice spacing, set by (3, and for bare 
quark masses appropriate to the b. Errors in the perturbative shifts are estimates of 
unknown higher order corrections. (Davies (1997)). 



from the perturbative calculation. It is also true that B for mqa — > oo vanishes 
at 0(a s ) so Z m — > 1. Since the dispersion relation for mesons at finite momen- 
tum cannot be obtained from potential model approaches, the quark mass there 
has to be fixed in a different way to the direct NRQCD method above. If the 
energies of states are calculated using the lattice potential including V c , rriQa 
should be adjusted until experiment is matched for, say, the T on applying a 
shift 2m Q a - V c (Bali, Schilling and Wachter (1997a)). 

Figures 15 and 16 show recent results for the bottomonium spectrum from 
lattice NRQCD. The errors shown on the plot are statistical errors only - it is 
clear that they are significantly smaller than those from light hadron calcula- 
tions. The simple form of the evolution equation for calculating the heavy quark 
propagator means that an average over a very large ensemble of gluon fields can 
be obtained with moderate computing cost. Also several different starting points 
can be used on a single gluon field configuration. 

The sources of systematic error are also under better control than for light 
hadron calculations. There, one of the most serious problems is that of finite 
volume. A large enough lattice is required not to squeeze the mesons under 
study and distort their masses. Because the Y is much smaller than, say, the 
p, a smaller space-time box is sufficient for its study. The calculations shown 
were done on lattices of around 1.5fm on a side. However, radial and orbital 
excitations are larger than ground states and such a lattice may be too small for 
35 and 2P states. Studies on larger volumes should be done for these in future. 
In fact direct calculations of the spectrum have worse finite volume errors than 
calculations of the heavy quark potential, because of lattice symmetries that 
protect V{R) (Huntley and Michael (1986)). 

Discretisation errors are an additional source of systematic error and in all the 
results shown, the leading discretisation corrections have been made as described 
above. 

Since NRQCD is a non-relativistic expansion, there are systematic errors 
from higher order relativistic terms that have been neglected. For the spin- 
independent terms all groups have used the lattice-discretised version of the 
Lagrangian of equation 22. This includes leading terms of 0{m,QV 2 ) and correc- 
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Fig. 15. Radial and orbital excitation energies in the X spectrum from lattice NRQCD 
obtained by different groups (Davies et al (1997), SESAM (1997)). Errors shown are 
statistical only. The scale has been set by an average of a -1 from the 2S — IS splitting 
and that from the IP — IS. Experimental results are given by dashed lines and for the 
1 Pi state is taken as the spin-average of the Xb states. 

tions of 0(rriQV 4 ). Errors are therefore at 0(niQV 6 ), giving v 2 x (a typical kinetic 
energy) = 0.01 x 400 McV for bottomoniuni = ~ 4 MeV. This error would be 
invisible on Figure 15, since it is a 1% error in the splittings shown. For the spin- 
dependent terms a 4 MeV error is more significant since splittings are smaller. 
This is the error estimate if only the leading spin-dependent terms of equation 
22 are used, as by the NRQCD collaboration (Davies et al (1994a)). The SESAM 
group (results presented at this school by Achim Spitz - see Spitz (1997) and 
SESAM (1997)) however has used the additional relativistic corrections to the 
spin-dependent terms given in the continuum by (Lepage et al (1992)): 



SH = -c 7 



9 



{dV-b} 



(60) 
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Cg 



35 



64m c 



■{D 2 ,cr • (DxE-ExD)} 



(T-ExE. 



This could reduce the relativistic error by another factor of v 2 to around 1% for 
spin splittings also. The SESAM group tadpole- improve all the terms above using 
uql and set the Cj to 1. In principal, however, unknown radiative corrections from 
lower order terms, e.g. 0(a s ) corrections to C4 beyond tadpole-improvement, can 
produce errors at the same order as the terms of (60) (if we take a s <~ v 2 ~ 0.1, 
but see Bodwin et al (1995)) so this is not a complete calculation at the next 
order. Note that relativistic corrections to spin-dependent terms are not known 
for a potential model. 



MeV 
20-- 
0-- 
-20-- 
-40-- 



hb 



ft 



Vb 



. fit" Xw 
«-f Xbi 



itt 



XbO 



Experiment 

Manke et al (UKQCD) (n f = 0),/3 = 6.0, v 6 , u 0P 
NRQCD (n f =0),/3 = 6.0, v 4 , u 0P 
NRQCD (nf =2, iTS, am, = 0.01), /3 = 5.6, v 4 , u 0P 
SESAM (nf = 2, W, n = 0.157), (3 = 5.6, w 6 , u 0L 
SESAM (nf = 2, W, n = 0.1575), /3 = 5.6, v G , u 0L 



Fig. 16. Fine structure splittings in the ground state T spectrum from lattice NRQCD 
obtained by different groups (Manke et al (1997), Davies et al (1997), SESAM (1997)). 
Errors shown are statistical only. The scale is set as in Figure 15. Experimental results 
are given by dashed lines. The spin-average of the Xb states has been set to zero. 



Figure 15 compares radial and orbital splittings to experiment. The lattice 
spacing chosen to set the scale is an average of that from the IP — IS splitting 
and that from the 2S — IS splitting. One striking feature is the disagreement 
with experiment for the calculation on quenched configurations (n/ = 0). The 
results on the partially unquenched configurations (n/ = 2) give much better 
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agreement. A quantity which exposes this error in the quenched approximation 
is the ratio of the 25*— 15 splitting to that of the IP— IS 1 . Figure 17 demonstrates 
that the fact that this ratio is too large is a physical effect - it is not affected 
substantially by lattice discretisation errors. We expect such an effect in the 
quenched approximation because a s runs incorrectly between scales. This means, 
as discussed earlier, that the quenched heavy quark potential is too high at small 
values of R and so the 5 states are pushed up with respect to P states, making 
the IP — 15 splitting too small relative to the 25— 15. The effect may be stronger 
for higher excitations, but they are subject to larger lattice errors. 

The error from the quenched approximation is even bigger if we look at quan- 
tities sensitive to a larger disparity of scales to maximise the effect of incorrect 
running. Figure 18 shows the ratio of the IP— 15 splitting in bottomonium from 
the NRQCD collaboration to the p mass from the UKQCD (UKQCD (1997)) 
and GF11 collaborations (Butler et al (1994), Weingarten (1997)). The UKQCD 
p mass results have included discretisation corrections for the light quarks; the 
GF11 results have not. It is clear that, although a result independent of lattice 
spacing is obtained when the p mass is improved, it is wrong. 
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Fig. 17. Dimensionless ratios of various splittings to the xii — T splitting against the 
lattice spacing in fm, set by the Xb — T splitting, in the quenched approximation. 
Experimental values are indicated by lines. (Davies (1997)). 



The fine structure in the spectrum is shown in Figure 16. As already dis- 
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Fig. 18. Dimensionless ratio of the yj, — T splitting to the p mass against the lattice 
spacing in fm, set by the Xt> — T splitting in the quenched approximation (Davies (1997)). 
Circles show the UKQCD improved p mass and the squares the GF11 unimproved mass. 
Experiment is shown by the burst. 



cussed, this is harder to calculate accurately than the spin-independent spectrum 
because it only appears as a relativistic correction. There is no large leading or- 
der term with non-perturbatively determined coefficient to stabilise the results 
as there is in the spin-independent case (D 2 /2toq). Since the clover discretisa- 
tion of E and B fields each contain four links (see figure 10) there are several 
powers of u in each term when tadpole-improvement is undertaken. This means 
that spin-splittings are affected strongly by the value of u and if u were set to 
1 (i.e. no tadpole-improvement) results much smaller than experiment would be 
obtained (Davies et al (1994a)). This also means that the spin splittings change 
when different definitions of uo are used, in the absence of a perturbative calcu- 
lation of the remaining radiative corrections (for preliminary results on these, sec 
Trottier (1997b)). The difference between uop and uql results in 10-20% shifts to 
the splittings at these lattice spacings (SESAM (1997)). The hyperfinc splitting 
is particularly sensitive; in leading order perturbation theory it is proportional 
to c\ (equivalent to Uq 6 when w factors in the D 2 term are taken into account). 
The presence or absence of the higher order relativistic corrections also affects the 
results at a similar level (Manke et al (1997), SESAM (1997)). Discretisation cor- 
rections, not surprisingly, are important as well. Because the fine structure (and 
particularly the hyperfine splitting) is sensitive to short-distance scales (consider 
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the functional form of the spin-dependent potentials), the discretisation errors 
can be quite severe. In the calculation of the NRQCD collaboration in which 
only leading order spin terms are included, the hyperfine splitting in MeV shows 
strong dependence on the lattice spacing - see Figure 19. This makes a physi- 
cal result hard to determine. Results of the other groups in Figure 16 include, 
along with the relativistic spin-dependent corrections, discretisation corrections 
to the leading spin-dependent terms (Lepage et al (1992)). This should reduce 
the lattice spacing dependence of the physical results but this analysis is not yet 
complete (see Figure 19 and Manke (1997)). Finally the spin splittings depend 
quite strongly on the quark mass (particularly again the hyperfine splitting) and 
for these the quark mass must be tuned accurately. This requires a very accurate 
determination of the meson kinetic mass as well as of the lattice spacing. 
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Fig. 19. The T — rib splitting in physical units vs the square of the lattice spacing in fm 
in the quenched approximation. The 2S — IS splitting has been used to set the scale. 
Squares are results from the lowest order action (Davies (1997)) and diamonds from an 
action which includes spin-dependent relativistic and discretisation corrections (Manke 
et al (1997), Manke (1997)). 



To compare to the real world we would like results with dynamical fermions 
of appropriate number and mass. To estimate how many dynamical fermions are 
'seen' by the bottomonium system, we need to know what the typical momenta 
being exchanged by the heavy quarks are. For T this qr is about 1 GeV, not 
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enough to make a cc pair, so the effective value of n/ should be 3. Almost all 
available sets of gluon field configurations have n/ set to 2, so extrapolation is 
necessary. The results will also depend on the light quark mass of the dynamical 
flavours, and this dependence at worst should be linear (Grinstcin and Rothstein 
(1996)): 



AM ~ AM 



1 + cV^. 

u.d,s 



(61) 



For qr ^> m q the answer with u, d and s quarks can be reproduced by 3 de- 
generate dynamical quark flavours of mass m s /3 (ignoring m u and mj). Results 
should be extrapolated as a function of dynamical fermion mass to m s /3 and 
then extrapolated to n/ =3 from n/ = and 2. 

In fact no significant ttlq dependence has been seen by the two groups, 
NRQCD and SESAM, who have done calculations on dynamical configurations 
(using different lattice formulations of the light fermions) . The SESAM collabo- 
ration with 3 dynamical quark masses does see a definite trend in m q , however 
(SESAM (1997)). Figure 20 shows the dependence of the ratio of the 25-15 
to IP — 15 splitting on n/ for the two groups. The results are consistent with 
experiment for n/ = 3 but rif — 2 cannot be ruled out without better statistics. 
More points at other values of n/ would be useful. 
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Fig. 20. The ratio of the T — T splitting to the \b — Y splitting as a function of the 
number of dynamical flavours. Circles are the results from the NRQCD collaboration 
(Davies et al (1997)) and the burst from the SESAM collaboration (SESAM (1997)). 
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The n/ extrapolation of the fine structure will be more difficult, even once 
physical results at a given n / are obtained. Since the fine structure probes much 
shorter distances it is possible that the effective number of flavours that it 'sees' 
is higher. Then the challenge will be to find appropriate quantities to set the 
scale for an extrapolation to, say, rif = 4. It will not be possible to use spin- 
independent splittings for which the real world rif value is 3 (Davies (1997)). 

Figure 16 shows disagreement with experiment for the P fine structure in 
the quenched approximation, both in overall scale and for the ratio p, equation 
37. Agreement is better on unquenched configurations, but the systematic errors 
described above must be removed before this is clear. The hyperfine splitting is 
very sensitive to the presence of dynamical fermions. It increases by <~ 30% as n f 
is increased from to 2 for the NRQCD results (see Figure 16). Extrapolations 
in rif using a variety of other short-distance quantities to set the scale (so the 
physical results differ from those in Figure 16) give a 'real world' value for the 
hyperfine splitting of around 40 MeV. The error is very large at present (25%) 
because of the inaccuracies in the fine structure. With improved calculations this 
can be reduced to 10%. 

2.4 Direct measurement of the charmonium spectrum on the lattice 

Unfortunately the NRQCD programme as described for bottomonium does not 
work as well for charmonium. It has been clear all along that charmonium is much 
more relativistic; with the NRQCD approach we can directly see the effects of 
higher order relativistic corrections to the Lagrangian. A calculation with the 
Lagrangian of equation 22 has errors at 0(rriQV 6 ) as discussed earlier (Davies et 
al (1995b)). This gives an error of around 30MeV which is 30% for spin splittings. 
On adding higher order terms these large corrections to the fine structure become 
manifest (Trottier (1997a)) and are actually rather worse than the naive 30%. 
An accurate calculation of the ^ — r\ c splitting, for example, would then require a 
high order in the NRQCD expansion, coupled with the determination of radiative 
corrections to the coefficients. 

It seems more useful to treat the c quark as a light quark and use stan- 
dard lattice approaches for relativistic quarks (Weingarten (1997), Montvay and 
Minister (1994)). However, the fact that m c a ~ 1 on current lattices can lead 
to significant discretisation errors. The heavy Wilson approach (El-Khadra et 
al (1997)) is an adaption of the standard Wilson light fermion action in which 
higher dimension operators are added to better match to continuum QCD by 
reducing errors of the form {pQa) n . The coefficients of these operators must be 
calculated and in the strict heavy Wilson approach they are considered as a 
perturbative series in a s but to all orders in niQa. At large niQa the Lagrangian 
becomes NRQCD-like (since no symmetry between space and time directions is 
imposed) and at small mqa it reduces to the form used for the Symanzik im- 
provement of light quarks. In principle it can span the region from one extreme 
to the other. In practise, NRQCD is rather simpler and faster to implement for 
really non-relativistic fermions. 
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The lowest order heavy Wilson action is identical to the Shcikosleslami- 
Wohlert (SW) action for light quarks in which a clover term 

Icffi^a^fXx) (62) 

with mga-independent coefficient is added to the Wilson action (Shcikholeslami 
and Wohlcrt (1985), Heatlie et al (1991)): 

S = J2^^n -K-Y^bPA 1 - l^Un^n+H + V'n+pl 1 + l^n^n]- (63) 
n n,/i 

The coefficient of the clover term is usually taken as 1 after the gauge fields 
have been tadpole-improved, and perturbative calculations to 0(a s ) (rriQa — ► 0) 
indicate no large additional radiative corrections (see Liischer et al (1997) for 
a discussion). Non-perturbative determinations of this clover coefficient for the 
light quark case are described by Sommer (1997). 

The meson energy-momentum relation must also be considered carefully for 
charm systems (El-Khadra et al (1997)). For the SW action there is an energy 
shift between the energy at zero momentum (the pole mass) and the kinetic 
meson mass that sits in the denominator of the kinetic term, as in the NRQCD 
case (equation 58). The shift increases as the quark mass increases and for niQa > 
0.5 it is important that the physical meson mass is taken from the kinetic mass 
and not from the pole mass which is used for light hadrons. ft is possible to 
remove this shift by adjusting coefficients in the full heavy Wilson approach, 
but it is not necessary. 

For the SW action there is also a problem with non- universality of the shift. 
It should appear simply as a shift per quark and therefore twice as big for a 
meson with two heavy quarks as for a meson with one heavy and one light 
quark. However there is a discrepancy between the shift per heavy quark in 
these two cases and it increases significantly as the heavy quark mass increases 
(Collins et al (1996a), Aoki et al (1997), see Figure 21). The discrepancy arises 
from lattice discretisation errors in relativistic D 4 -type terms in the heavy quark 
action affecting the heavy-heavy mesons. These terms need to be correct in order 
for the binding energy of a meson to be fed into its kinetic mass. Since the kinetic 
mass appears at 0(niQV 2 ) it is 0(mQV 4 ) terms which do this, whereas the pole 
mass is 0(1). If the C(toqw 4 ) terms are incorrect, the binding energy will appear 
only in the pole mass and the shift will then depend on the binding energy. For 
a heavy-light meson the binding energy is provided by the light quark and this 
problem does not arise. In the NRQCD case (Davies et al (1994b)) the relativistic 
terms are correct because the D 4 /8toq term is added by hand and discretisation 
corrections remove Df rotational non-invariance. For the SW action this isn't 
true; the D 4 term has an incorrect mass and there is an uncancelled Df term 
(Kronfcld (1997)). Both these effects can be corrected for in the heavy Wilson 
approach (El-Khadra et al (1997)), but this has not been done as yet. For charm 
quarks with m c a < 0.5 there is not a significant problem, as in clear from Figure 
21. 
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Fig. 21. The energy shift per heavy quark for the SW action, as a function of heavy 
quark mass, m^a. Note the difference between results for heavy-heavy mesons (open 
symbols) and those for heavy-light mesons (filled symbols). (Aoki et al (1997)). 



The calculation of charm quark propagators with the SW action proceeds as 
for conventional light quarks (Weingarten (1997), Montvay and Minister (1994)) 
with the calculation of rows of the fermion matrix inverse by an iterative pro- 
cedure. This converges quite rapidly for heavy quarks, but care must be taken 
to allow enough iterations for the solution to propagate over the whole lattice. 
Meson correlation functions are put together using various smearing functions 
and then fitted to multi-exponential forms as described for the NRQCD case. 
The charm quark mass is tuned by the kinetic mass method as before. 

Figure 22 shows the spectrum for charmonium calculated recently on quenched 
gluon fields with this method and presented at this school by Peter Boyle (Boyle 
(1997b)). Previous results (El-Khadra and Mertens (1995)) are in agreement 
with this, but don't show such complete fine structure. The hyperfine splitting 
is clearly underestimated and this could be a quenching error, since this splitting 
increases with n/, as discussed in the bottomonium case, or it could mean that 
c s is underestimated. The fine structure also shows a discrepancy for the ratio p 
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(see equation 37). 

The systematic errors of the SW action for charmonium need some analysis 
(El-Khadra and Mertens (1995)) before we can extract physical (quenched) re- 
sults from these calculations. Only one calculation on unquenched configurations 
has been done (Collins et al (1996a)) and problems with fitting errors made it 
hard to draw conclusions. 



GeV 



4- 



3.0 



Vc J/ijj XcO Xcl Xc2 h c 

Fig. 22. The charmonium spectrum from quenched lattice QCD using the SW action 
(Boyle (1997b)). 



Future work will need to investigate the use of actions with higher order 
terms, possibly on anisotropic lattices (Lepage (1996)). Some preliminary work 
on the charmonium spectrum has been done with these (Alford et al (1997)). 
A small lattice spacing in the time direction is useful for improving exponential 
fits, particularly for excited states, and does not need to mean a small lattice 
spacing in the spatial directions. Indeed such an anisotropy is very natural for 
non-relativistic systems, as we have seen. 



2.5 Other heavy-heavy states 

There is a lot of interest in the literature in other heavy-heavy bound states, 
which have not yet been seen experimentally. The one most likely to be seen in 
the near future is the mixed bound state of bottom and charm quarks; indeed 
candidates for the 1 5o ground state, the B c , have been seen recently (DELPHI 
(1997a), ALEPH (1997)). 

The be system actually has a lot in common with the heavy-light systems of 
the next section, although it is classified here as heavy-heavy because of its quark 



44 Christine Davies 



content. The charm quark in the B c will be more tightly bound and therefore 
more relativistic than in charmonium, and we have already seen that a non- 
relativistic approach to charmonium is rather inaccurate. In addition, because 
charge conjugation is not a good quantum number, the two 1 + P states will mix 
to give a different P fine structure to that for heavyonium. 

Recent continuum potential model results for the be spectrum are given in 
Eichten and Quigg (1994) and Gershtein et al (1995). 2 sets of S states, 1 set 
of P states and 1 set of D states are expected below threshold for the Zweig 
allowed decay to B,D (7.14 GeV). Note that the be states below threshold are 
particularly stable since the annihilation mode to gluons is also forbidden. 

First lattice calculations (Davies et al (1996)) have used NRQCD for both 
the c and b quarks. Agreement with potential model results was found within 
sizeable systematic uncertainties. Better calculations will use NRQCD for the b 
quark and relativistic formulations for the c quark (Shanahan (1997)). However, 
uncertainties still remain about how to fix the bare quark masses in the quenched 
approximation, because of the variations possible in the determination of the 
scale. These problems should become more tractable when complete calculations 
are done including the effects of dynamical fermions. Only preliminary results 
are available on unquenched configurations using NRQCD for b and c (Gorbahn 
et al (1997)). Figure 23 shows a comparison of the spectrum for lattice and 
potential model results. 

Other states of a more speculative nature are hybrid states which include a 
gluonic valence component; QQg. Observation of these states would be a direct 
confirmation of the non-Abelian nature of QCD (see possibly E852 (1997)). 
Hybrids are expected containing all quark flavours, but the advantage of studying 
the heavy-heavy hybrids is that the normal QQ spectrum is relatively clean, both 
experimentally and, as discussed here, theoretically. 

There have been several lattice calculations of the hybrid potentials for a po- 
tential model analysis of these states (see, for example, Perantonis and Michael 
(1990), Juge, Kuti and Morningstar (1997)). Wilson loops are calculated whose 
spatial ends have the appropriate symmetries to project onto the different hybrid 
sectors (see Figure 24). The hybrid potentials obtained can be compared to ex- 
pectations from excited string models and from bag models. From a Schrodinger 
equation, masses for the hybrid states can be determined; the particular interest 
is in 'exotic' states which cannot appear in the usual QQ sector. These states 
have quantum numbers (J = odd) h , (J = even) + ~ or 0~~. They are most 
likely to be visible if their energies are below threshold for Zweig-allowcd decay, 
but it is not clear where this threshold is. Some models expect the hybrid states 
to decay to an S'-wave heavy-light state and a P-wave, in which case the thresh- 
old is rather higher than for conventional heavyonium decay (for a recent review 
of expected hybrid phenomenology see Close (1997)). 

The hybrid potentials obtained (see Figure 25) are very flat, indicating broad 
states, closely packed in energy. The lightest mass hybrids from these potentials 
are close to the threshold described above. The same picture is obtained by 
calculating the masses of heavy-heavy hybrids directly using NRQCD (Collins 
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Fig. 23. Lattice results for the spectrum of the B c system. Open circles indicate 
NRQCD results on quenched configurations, closed circles those on partially un- 
quenched configurations from the MILC collaboration. Bowties indicate results on 
quenched configurations using NRQCD for the b quark and relativistic c quarks. Er- 
ror bars are shown where visible and only indicate statistical uncertainties. (Davies et 
al (1996), Shanahan (1997), Gorbahn et al (1997)). Dashed lines show results from a 
recent potential model calculation (Eichten and Quigg (1994)). 



Fig. 24. An example of an operator used at the end of a Wilson loop in the calculation 
of the hybrid (77) potential. 



(1997b), Manke (1997)). Further work must be done on the spectrum if the 
states are to be accurately predicted for experimental searches. 
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R/r„ 

Fig. 25. The heavy hybrid potential from a recent lattice calculation (Juge, Kuti and 
Morningstar (1997)). The Eg potential is the usual central potential; the IJ potential 
has a gluonic excitation with spin 1 about the QQ axis and the A potentials have spin 
2. 

3 The heavy-light spectrum 
3.1 Mesons 

These are bound states with one heavy valence quark or anti-quark and 1 light 
anti-quark or quark. The levels show a similar picture to that for the heavy- 
heavy spectrum with the lightest state the pseudoscalar ( 1 5o) and close by the 
vector ( 3 Si). For charm- light we have pseudoscalars cd = D + , cu = D° and cs 
= D s , and vectors, D*°, D* + and D* s . For bottom-light we have pseudoscalars 

bd = B , bu = B~ and bs = B s , and vectors again for each. We shall ignore the 
distinction (and slight mass difference) between heavy-light mesons containing 
u and d quarks and often just refer to D and B. Radially excited S states, D' 
and B 1 are about 500 MeV above the ground states and below these come a set 
of positive parity P states, denoted by their spins Dq/Bq, Dx/Bi, D^fB^ etc. 
or more generically D** and B**. See Figure 26 and The Particle Data Group 



(1996). 
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Fig. 26. The spectrum of heavy-light mesons. 



We do not expect a potential model to work well for the heavy-light spectrum 
because the light quarks arc now relativistic. The heavy quarks are still non- 
relativistic, however. Taking Aqcd as a typical QCD momentum scale of a few 
hundred MeV, we have 

MomentuniQ <~ Momentum^ ~ Aqcd 
vq Aqcd 
c rriQ 

giving vq ~ 0.1 for b in B and 0.3 for c in D. This is vq, not Vq, so the heavy 
quark is actually more non-relativistic than in heavy-heavy systems. 

A useful analysis is provided by Heavy Quark Symmetry (see Neubert (1994) 
for a review). In the itiq — > oo limit QCD has an [SU(2Nf)] symmetry where 
Nf is the number of flavours of heavy quark in the theory. This is evident 
from the NRQCD Lagrangian of equation 22 which, by the arguments above, is 
appropriate for the heavy quarks here also. As rriQ — > oo the Lagrangian becomes 
^ D t 4> when the quark mass term is removed. Thus the heavy quarks become 
spinless and any distinction between flavours disappears (apart from the overall 
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energy level set by the missing mass term). The picture of a heavy- light meson 
becomes one of a static heavy quark surrounded by a fuzzy cloud of the light 
degrees of freedom, known as 'brown muck', as in Figure 27. Interactions that 
probe only momentum scales appropriate to the brown muck will not be able 
to see details of the heavy quark at the centre. Notice that this is a completely 
different physical picture to that for heavy-heavy mesons. 




Fig. 27. A heavy-light meson in the Heavy Quark Symmetry picture. 



From the picture of Figure 27 there is a natural distinction between energy 
shifts in the spectrum that are caused by something changing for the light degrees 
of freedom, e.g. a radial or orbital excitation, and those caused by something 
changing for the heavy quark, such as its flavour or spin. In the first case we ex- 
pect that radial and orbital excitation energies should be approximately indepen- 
dent of the heavy quark flavour, and in the second case we expect much smaller 
splittings with strong heavy quark mass dependence between states of different 
Sq. This hierarchy of splittings is similar to that for heavy- heavy mesons but for 
different reasons. The mQ-indcpendcncc of the IP — IS splitting in heavyonium 
is an accident; in heavy- light mesons it is the consequence of a symmetry of the 
non-relativistic effective theory as uiq — > oo. 

A power-counting analysis of the terms in the NRQCD Lagrangian is useful 
to demonstrate this effect (Ali Khan et al (1996)). 



D t ~ D ~ Aqcd 



(64) 
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This shows that, for the heavy-light case, the NRQCD Lagrangian is a 1/mg 
expansion (unlike the heavy- heavy case where terms at different order in 1/ uiq 
appeared at the same order in Vq). The leading order term is the D t term and 
then at the next order come two 1/mg terms - the kinetic energy of the heavy 
quark and the spin coupling to the chromo-magnetic field. These are the first two 
terms to know about the heavy quark flavour (mass) and its spin. Any splitting 
that requires this knowledge will appear first at 1/rag in an expansion in the 
inverse heavy quark mass. 

The heavy quark spin, Sg, is a good quantum number in the heavy quark 
limit and so we can classify states according to j; — J — Sq. Each j; state becomes, 
on the addition of the heavy quark, a doublet with J = ji ± 1/2 (Isgur and Wise 
(1991)). An analogy can be drawn with atomic physics and the decoupling of 
the nuclear spin as m e /mjv — * 0. The lightest states are the L = 0, ji = 1/2, 
3 Si (D*,B*)/ 1 S (D, B) doublet. For heavy-light P states the light quark spin, 
S q , is coupled to the orbital angular momentum to make states of overall spin, 
ji = 1/2(2 polarisations ) or ji = 3/2 ( 4 polarisations). Coupling Sq to ji = 
3/2 gives total J=2 {B*,D*) and J=l (B[, D[) (8 states altogether). Coupling 
Sq to ji — 1/2 gives J=0 or J=l (4 states altogether). Thus in the jj coupled 
basis we reproduce the same 12 states as the LS coupled 1 Pi, 3 Po,i,2 multiplet. 
However, the spin 1 states are a mixture of the 1 P i and 3 Pi (with mixing angle 
35°) because of a lack of charge conjugation. In the toq — ► oo limit only the 
splittings caused by the light degrees of freedom remain. The jj basis becomes 
the correct one and all the ji = 3/2 states become degenerate but split from 
all the ji = 1/2 states. The ji = 3/2 states are narrow (and therefore visible) 
because of the high orbital angular momentum required in decays to D^B^n 
for J=2. J=l can only decay to D* /B*ir but, having the same ji as the J=2 
state, has a similar total width (see Figure 28 and Isgur and Wise (1991)). 

The difference between the ji + 1/2 and ji — 1/2 members of a doublet is 
a spin flip of the heavy quark. The leading term that gives rise to this in the 
NRQCD Lagrangian is the <tq • B term, yielding a splitting behaving as Ax (spin 
factors) x 1 / toq . A is an expectation value in the light quark degrees of freedom so 
the heavy quark mass dependence of the splitting is as 1/mg in leading order. 
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Table 4 shows experimental values for the vector-pseudoscalar splitting (The 
Particle Data Group (1996)). 1/rriQ behaviour fits very well if we take m c ss 1.5 
GeV and mi, m 5 GeV. We can also consider the strange quark as a heavy quark 
rather than a light one and add the value for the strange- up/down system, the 
K into the Table. This only works moderately well, with m s « 0.5 GeV, say. The 
coefficient of the I/tuq dependence is of order 0.2 GeV 2 , which is compatible 
with an expectation value in a light system of Aq CD . Note that the variation 
with light quark mass between u/d and s is very small. It is clear that B c cannot 
be fitted into this heavy-light picture since its expected hyperfine splitting is 
much larger (see section 2.5). 



Splitting 


Experiment/ MeV 


'Expected' value / MeV 


K* — K 


398 


457 


D* - D 


141 


152 


D* s - D 3 


144 




B* -B 


46 


46 


Bt - B s 


47 




A" 2 *(1430) - #i(1270) 


154 






37 




D* s2 - D sl 


38 





Table 4. Hyperfine splittings for different heavy-light systems; the top group for the 
L = 0, ji = 1/2 doublet, the lower group for the L — l,ji — 3/2 doublet. The final 
column gives expected values for the 3 Si — 1 So splitting rescaling from the B system by 
the inverse ratio of quark masses given in the text. (The Particle Data Group (1996)). 



Table 4 also shows results for the L = 1, ji = 3/2 doublet from the D and 
K systems (Eichten et al (1993)). The ratio of splittings between D and K is 
rather different for this case to the one above, probably showing that the K 
is stretching the limits of HQS arguments. Nevertheless, we expect a splitting 
B* 2 -B x of m c /m h x (D^-D^ - 12 MeV A value of 26 MeV is given in DELPHI 
(1995) for the B s . Experimental results for the L = 1 ji = 1/2 doublet arc not 
available since these are much broader than the ji = 3/2 doublet. 

In contrast there are several splittings that we expect, from the arguments 
above, to be controlled by changes in the light quark degrees of freedom and 
therefore to be independent of the heavy quark mass at leading order. One of 
these is the splitting between the heavy-strange and heavy-up/down mesons. 
Experimentally this is satisfied at the 10% level (see Table 5). Other such split- 
tings are those between radially excited S states and the ground states for which 
experimental information is very limited (B — B = 580McV (Landua (1997)) 
and D*' - D* = 630 MeV (DELPHI (1997b))), and between orbitally excited P 
states and the ground 5* states, which we discuss below. 
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Splitting 


Experiment / MeV 


D s -D 


99 


B 3 ~B 


90 



Table 5. Experimental values for splittings between heavy-strange and heavy-up/down 
systems (The Particle Data Group (1996)). 

\ — J=2 

J=l l/mq 

a Pi , 

'Spj' nriQ — ► oo 



3 P J=l } 1 / m Q 

31 = 1/2 



Fig. 28. On the left, P fine structure for a degenerate heavy-heavy system. On the 
right, P fine structure for a heavy- light system. 

For the orbital splittings between heavy-light S and P states we should cal- 
culate a splitting between spin-averaged states, to remove the spin-dependent 
1/toq effects, and make as clear as possible the mg-independent light quark 
effects. Since jt = 1/2 states have not been seen, we compare in Table 6 the 
splitting between the spin- average of the ji — 3/2 P states and the spin- averaged 
S states for D and B. Good agreement between the 2 systems is seen. For B 
states good spin separation of the P states is not yet available. 

We would also expect the splitting between the ji = 1/2 and the ji = 3/2 
states to be approximately independent of ttiq (Isgur (1997)), although the phys- 
ical spin 1 states will be a mixture of the jj states away from the toq — ► oo limit. 
This cannot be checked experimentally as yet. 

3.2 Baryons 

There is a huge array of baryon states with one heavy quark and two light 
quarks. Again we can make sense of their masses using Heavy Quark Symmetry 
arguments. We view the baryon as a static colour source (for the heavy quark) 
surrounded by a fuzzy light quark system which is made of two light quarks this 
time instead of a light anti-quark (Falk (1997)). 



52 Christine Davies 



Splitting 

if p (1368) -¥(792) 
~5 P (2445) -~D{ 1975) 
i) sp (2559) -i3 s (2076) 
B**(5698) -5(5313) 
B;*(5853) - B s (5404) 



Experiment / MeV 
576 
470 
483 
385 
449 



Table 6. Splittings between spin- averaged ji = 3/2 P states (or experimentally un- 
separated P states) and spin-averaged S states for heavy-light systems (The Particle 
Data Group (1997)). 



We will discuss only the case of zero relative orbital angular momentum. For 
two different light quarks we can combine the light quark spins to give a total Si 
of or 1 . If the light quarks have the same flavour, only Si = 1 is possible by Fermi 
statistics, remembering the overall anti-symmetry of the colour wavefunction. 
Coupling the heavy quark spin then gives the combinations in Table 3.2, with 
overall spin-parity assignments. 



baryon Qqq Si J p mass c /MeV masst/MeV 



A Q[ud] 0+ Y~ 2285(1) 5624(9) 

E Q{ud},uu,dd 1+ i + 2453(1) 5797(8) 

E* Q{ud},uu,dd 1+ | + 2519(2) 5853(8) 

S Q[u/ds] 0+ | + 2468(2) 

5' Q{u/ds} 1+ i + 2568(?) 

S* Q{u/ds} 1+ | + 2645(2) 

Q Qss 1+ i + 2704(4) 

3 + 



Q* Qss V 



2 



Table 7. J p possibilities for baryons containing one heavy quark along with two light 
quarks. The names are given with subscripts c or b. Masses are given in the last two 
columns, taken from The Particle Data Group (1997), DELPHI (1995) for E b and 
WA89 (1995) for s' . 



Using HQS arguments we would expect the splitting between the spin average 
of S and S* states and the A to be independent of rriQ, since this splitting 
represents a change in ji. We can check this in Table 8, and it works well even 
when the s quark is considered as a heavy quark. There is in fact very little room 
for sub-leading 1/Mq dependence which can in principle be there (A^ CD /m c <~ 
50 MeV) . In the last row is given for comparison the splitting between the spin- 
average of the E* and E c and the E c . This is essentially the same splitting except 



The Heavy Hadron Spectrum 



53 



for the different light quark content. The answer is significantly different, showing 
more sensitivity to light quark content than for the mesons (Falk (1997)). The 
physical E c and E c will be mixtures of the HQS states, just like the spin 1 
meson P states, but this should not be a big effect. An equal spacing rule, 
Q c — E c = E c — S c holds well. 



Splitting 


Experiment /MeV 


£~ a - A a 


203 


S c -A c 


212 




210 


„ 

" c 


150 



Table 8. Experimental values for the splitting between the spin-average of S states 
and the A for different heavy quarks including s. In the last row a comparable splitting 
is given for the S c . 



All the fine structure splittings between states of the same Si but different J 
should behave as 1/mg. Table 9 shows the experimental information on this for 
the £ baryons. The s quark fits well into this picture, but the experimental Sf } — 
E splitting looks significantly different from the expected value (Falk (1997)). 
The experimental results need to be confirmed, however. The E* — E splitting 
agrees well with the S* — £ showing no large m s effects here. 



Splitting 


Experiment / MeV 


'Expected' / MeV 


S s — S s 


191 


198 


z* - s c 


66 


66 


£b — £b 


56 


20 




80 





Table 9. Splittings between S* and S states for different heavy quarks, including s. 
The last column gives expected values for I/tuq behaviour compared to the splitting 
for c. 



We can also take splittings between baryons and mesons. The simplest split- 
ting is between the A baryons and the S state mesons. To remove spurious toq 
dependence we should take the spin- average of the 1 So and 3 5i meson states 
(Martin and Richard (1987)). Table 10 shows the experimental results; again 
Heavy Quark Symmetry works much better than might be expected. 

HQS yields only the vhq dependence of the splittings; it must be combined 
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Splitting 
A s -K_ 

A c -D 
A b -B 



Experiment / MeV 



323 
310 
310 



Table 10. The splitting between the A baryon and the spin average of S state 
heavy- light mesons for different heavy quarks, including s. 



with a non-perturbative method of determining the coefficients of this depen- 
dence. QCD sum rules can be invoked here (Neubert (1994)); Lattice QCD 
provides a better ab initio method. We discuss results from lattice QCD in the 
next subsection. 

Exercise: Discuss what you would expect for heavy-heavy-light baryons. 
Take Qi ± Q 2 . 

Exercise: Compare orbitally excited A s and A c baryons from the Par- 
ticle Data Tables. What does this lead you to expect for the orbitally 
excited A b ? (Rosncr (1995)). 

3.3 Direct calculations of the heavy-light spectrum on the lattice 

Following the methods described for the heavyonium spectrum, we can calculate 
the heavy-light spectrum directly using lattice QCD. We must combine a heavy 
quark propagator with a light anti-quark propagator (or two light quark propa- 
gators) to make a meson (baryon) correlation function. This we fit as before to 
a sum of exponentials to extract ground and excited state energies and masses. 
In principle some of the excited states can undergo strong decays upsetting this 
relation, but this does not happen in current lattice simulations. 

For hadrons containing a b quark the best method is probably to use NRQCD 
for the heavy quark as described for bottomonium in section 2.3. Because of the 
different power-counting rules for the heavy-light case, the Lagrangian used can 
be different to that for heavyonium. For example, a consistent calculation to 
0(l/mg) would include D t , T> 2 /2niQ and cr ■ H/2rriQ terms (tadpole-improved 
as before). In fact for the heavy- light case the spectrum can be calculated in 
the static limit with simply the D t term, because the light quark provides the 
kinetic energy. In this case, of course, only states of a given ji are obtained 
with no hypcrfine splittings. The static limit is very cheap computationally but 
much noisier (Lepage (1992)) than NRQCD even at very large mg and for this 
reason it may be more accurate to obtain static results from the limit of NRQCD 
calculations. For hadrons containing a c quark, we will discuss results using the 
heavy Wilson (SW) action. 

Since we do not have a potential model in principle to guide our intuition, 
it is more difficult to think of good smearing functions for heavy-light mesons. 
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Fig. 29. The B spectrum from lattice QCD using NRQCD for the b quark. Circles 
are in the quenched approximation; open circles use m s from K and closed circles, 
ms from K* . Squares are results on configurations with n/ = 2 dynamical fermions. 
Experimental results (The Particle Data Group (1997)) are given by dashed horizontal 
lines. The B meson mass is fixed to its experimental value in all cases (Ali Khan 
(1997)). 



A lot of effort has been put into this for mesons in the static case (Duncan 
et al (1995), Draper et al (1995)) to ameliorate the noise problems. Again, the 
smearing does not affect the values of masses obtained, but a good smearing can 
reduce the errors. In fact potential-model type wavefunctions (much broader 
than for hcavyonium) do work reasonably well (Duncan et al (1995), Ali Khan 
et al (1996)), used as a source for the heavy quark, as do gauge-invariant smear- 
ings typical of light hadron calculations (UKQCD (1996b)). The light anti-quark 
for the meson is taken to have a delta function source. Alternatively both prop- 
agators can be smeared, and for baryons it is certainly a good idea for all the 
propagators to be smeared (UKQCD (1996b)). 

The meson operators are similar to those for heavyonium - tpgil^Xq- For the 
NRQCD heavy quark case Q is a 2 x 2 matrix in spin space and only 2 components 
are taken from the 4-component light quark. The colors of heavy quark and 
light anti-quark are matched for a colour singlet. The baryon operators need an 
anti-symmetric colour combination, and the light quark propagators combined 
with appropriate spins (UKQCD (1996b)). For example the Aq operator (with 
smearing factors, (j> suppressed) is: 

O = e ABC ( X £ C 75X f 2 )Vj? ( 68 ) 

where C is the charge conjugation matrix. 

In principle, having calculated the bottomonium spectrum in NRQCD as 
in section 2.3 on a given set of gluon configurations, we can determine a -1 
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and the bare b quark mass, nib, and the calculation of the B spectrum should 
have no parameters to tune. Unfortunately this is not true in the quenched 
approximation. The disagreement with experiment shown in Figure 18 makes 
it clear that the a -1 fixed from the T spectrum would be ~ 20% different to 
that from M p , because of the different momentum scales appropriate to the two 
systems. Heavy-light systems are much closer to light hadrons in these terms than 
to heavyonium. For the best quenched results we really need to use a value for 
a -1 from the heavy- light system itself, but the lack of experimental information 
on P states makes this hard, since the obvious quantity to use is the IP — IS 
splitting. Usually a -1 is taken instead from light hadron spectroscopy. Large 
statistical and systematic uncertainties there then give a rather large error. 
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Fig. 30. Masses of baryons containing one b quark from lattice QCD. Circles use 
NRQCD for the b quark in the quenched approximation, the box uses NRQCD on 
configurations with n/ = 2 flavours of dynamical fermions. Triangles (Alexandrou et 
al (1994b)) use Wilson fermions, and diamonds the SW action (UKQCD (1996b)) ex- 
trapolating from the region of the charm quark, again in the quenched approximation. 
Experimental results are given by horizontal lines (The Particle Data Group (1997), 
DELPHI (1995)). (Ali Khan (1997)). 



This creates a problem with the bare b quark mass, m,b, since it was fixed 
in bottomonium using a -1 from that system. It should be fixed again in heavy- 
light systems using the kinetic mass of, say, the B. This is difficult to extract 
accurately because Bs are lighter than Ts. E{p) — E(0) is larger for the B 
than the Y so the noise in the meson correlation function at finite momentum 
p (set by £'(0)) is worse. An alternative is to calculate the usual energy at 
zero momentum, Eb(0), and apply the energy shift per quark in lattice units 
calculated for heavyonium to get mgo (Ali Khan et al (1996), Collins et al 
(1996b)). 
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These problems mean that the heavy-light spectrum cannot be as accurately 
calculated as that of heavyonium. Once dynamical fermions are included suffi- 
ciently well to mimic the real world there can only be one value of a -1 and m b / c . 
We are a long way from this point at present, however. It is not even possible to 
perform consistent rif extrapolations (to n/ = 3?) of the heavy- light spectrum 
from results at nj — and 2 (Collins et al (1996b) and in preparation). For 
heavyonium differences in methods of fixing a -1 disappeared on this extrapola- 
tion but this is not currently true for heavy-light mesons and shows the presence 
of systematic errors. 

Another difficulty with the heavy-light spectrum is that of fixing the light 
quark mass. This is a problem shared with light hadron calculations (Weingarten 
(1997), Montvay and Minister (1994)). The B and D calculations must be done 
with several different light quark masses far from the physical u/d masses and 
the results extrapolated to the chiral limit. This inevitably causes an increase 
in statistical and systematic errors. For the B s and D Sl it is possible to inter- 
polate to the s quark mass although there are ambiguities in fixing that, again 
possibly arising from the quenched approximation (see, for example, Gupta and 
Bhattacharya (1997)). 

Figure 29 shows the 6-light meson spectrum using NRQCD for the b quark, 
fixing nib from the B mass and a -1 from M p (from a recent review by Ali 
Khan (1997)). The overall agreement with experiment is good. The B* — B 
splitting is too small, however, both on quenched and on partially unquenched 
configurations. As before, this may be a quenching effect and/or it may arise 
from radiative corrections to C4 beyond tadpole-improvement (see the discussion 
for heavyonium). The problems with fixing m s are clear. The P states still 
have rather large error bars but the ordering, Bq , B\ , B\ is becoming clear in 
the lattice results (in disagreement with some expectations (Isgur (1997))). The 
spin 1 P states cannot be clearly separated as yet. Experimental results on the 
P states are likewise uncertain. Results at a different value of the lattice spacing 
are compared in Hein (1997). 

Figure 30 shows the 6-light baryon spectrum using NRQCD for the b quark 
(Ali Khan (1997)). Agreement with experiment is again reasonably good, al- 
though the Ab baryon is apparently too heavy on the partially unquenched con- 
figurations. The baryons are probably rather susceptible to finite volume effects, 
and further work is definitely needed on bigger volumes. The S* — £ splitting 
is too small in the quenched approximation, which does not seem surprising by 
now. 

Results in the static limit for mesons and baryons, for the states that still 
exist there, are similar to those from NRQCD but have been in the past less 
accurate - see Peisa and Michael (1997), UKQCD (1996a), Duncan et al (1995), 
Alexandrou et al (1994a) and Duncan et al (1993). A comparison is made in 
Figure 30 to results using heavy Wilson quarks for the b. This will be discussed 
further below. 

Arguments earlier showing that the heavy quarks are more non-relativistic 
in heavy-light than heavy-heavy does mean that NRQCD should work better for 
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Fig. 31. The D meson spectrum from lattice QCD using a tadpole-improved SW action 
in the quenched approximation. Masses are fixed relative to the spin average of the D a 
and the D* a (Boyle (1997a)). Horizontal dashed lines mark the experimental results 
(The Particle Data Group (1996)). 



the D than for the ip but currently results are only available for ^-states at one 
value of the lattice spacing (Hcin (1997)). For c-light mesons and baryons the 
SW action (or other heavy Wilson action) is probably to be preferred even in this 
case. Figure 31 shows a recent D spectrum using the tadpole-improved SW action 
for the c quark presented here by Peter Boyle (Boyle (1997a)). The agreement 
with experiment is encouraging but it has not been possible to extract all the P 
fine structure as yet, and errors bars there are still rather large. Uncertainties in 
how to fix a" 1 and m c are the same here as for the b case above. In the Figure 
a -1 is taken from light hadron spectroscopy. Results have been compared at two 
values of the lattice spacing (Boyle (1997b)). No D spectrum is available from 
unqucnchcd configurations as yet. 

Figure 32 shows the c-light baryon spectrum using the SW action for the c 
quark but this time not tadpole-improved (UKQCD (1996b)). Agreement with 
experiment for the mQ-independent splittings is reasonable but the hyperfinc 
splittings are much too small. At least a part of this comes from the lack of 
tadpole- improvement since this directly affects the '<tq • B' term in this formal- 
ism. 

It is tempting to try extrapolating to the b quark from the c quark using 
the results from the SW action and HQS arguments to set the mQ-dependence 
of splittings. This is probably fine for splittings which have little or no toq 
dependence, in which case extrapolation is not really necessary. We have seen 
that there are several splittings for which the leading behaviour is a constant 
and for which there seems almost no sub- leading dependence on vhq. For the 
splittings that have strong vhq dependence it is much more difficult to pick up 
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Fig. 32. The spectrum of baryons containing one c quark obtained on the lattice using 
the SW action for the c quark in the quenched approximation (UKQCD (1996b)). The 
horizontal dashed lines give experimental results (The Particle Data Group (1997)). 



this dependence from the small toq side than from the large. Figure 30 compares 
results from NRQCD with extrapolated results from unimproved Wilson quarks 
(Alexandrou et al (1994b)) and SW quarks (UKQCD (1996b)). In the latter 
case the low value obtained for the hyperfine S* — £ splitting becomes worse on 
extrapolation. In principle the SW action is safer for heavy-light mesons than 
for heavy-heavy as was discussed in section 2.4. and so calculations at the b itself 
can and should be done with this method (Simonc et al (1997)). 

Finally lattice QCD calculations do not have to restrict themselves to he 
physical quark masses but can explore the whole heavy quark region. This en- 
ables a fit to the dependence on vtiq (or on the pseudoscalar meson mass, say) of a 
range of splittings. The coefficients of this dependence are then non-perturbative 
parameters of a heavy quark expansion which can be made use of in other heavy 
quark relations. The values of the coefficients can be compared to expectations 
of powers of Aq C d an d to QCD sum rule results (Collins et al (1996b), Collins 
(1997a), Gimenez et al (1997)). 



4 Conclusions 

There has been a lot of progress in heavy hadron spectroscopy using the tech- 
niques of lattice QCD in recent years, converting the qualitative understanding 
of potential models and Heavy Quark Symmetry into clear numerical results 
that test QCD. 
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Further work is still needed to bring down systematic errors. Upsilon spec- 
troscopy is the most accurate at present. Here, more calculations need to be done 
with a non-relativistic action which includes next-to-leading spin-dependent terms 
and radiative corrections to leading terms. Finite volume effects must be studied 
for radially excited states. More accuracy is needed on dynamical configurations 
with several different values of n j to allow for a clear extrapolation of fine struc- 
ture to the real world. A prediction for the T — r/i, mass at the 10% level is a 
realisable goal with current calculations. 

The charmonium spectrum is more complete experimentally but more work 
is needed on the lattice using heavy Wilson actions to reduce statistical and 
systematic errors (e.g. from D 4 terms). Calculations on configurations with dy- 
namical fermions are required for extrapolations to compare to experiment. 

In the heavy- light sector statistical and systematic errors are inevitably larger 
and these must be reduced if we are to get a clear picture of the fine structure 
and radial excitations from the lattice that are now being seen experimentally. 
Analyses of scaling as the lattice spacing is changed and finite volume studies for 
these systems are still at an early stage. In the next few years a clearer picture 
will emerge of the effect of the quenched approximation on 'softer' momentum 
systems such as light and heavy-light hadrons and ambiguities of scale setting 
and quark mass fixing should be removed. 

Finally, as noted at the beginning, we are also interested in matrix elements 
for radiative and weak decays of heavy hadrons, particularly those which are 
important for the experimental B physics programme. Calculations of these are 
being done on the lattice also, using the techniques described here for the spec- 
trum. These calculations are much harder and accurate spectrum results will be 
a prerequisite for accurate matrix elements. 

Acknowledgements 

I thank the organisers for a very enjoyable school and the following for help 
in preparing these lectures: Arifa Ali Khan, Gunnar Bali, Peter Boyle, Sara 
Collins, Joachim Hcin, Henning Hoeber, Peter Lepage, Paul McCallum, Colin 
Morningstar, Junko Shigcmitsu, John Sloan and Achim Spitz. A lot of the work 
described here was supported by PPARC and NATO under grant CRG 941259. 
I am grateful to the Institute for Theoretical Physics, UCSB, for hospitality and 
to the Leverhulme Trust and the Fulbright Commission for funding while these 
lectures were being written up. 

References 

ALEPH collaboration (1997): Phys. Lett. B402, 213. 

C. Alexandrou, S. Giisken, F. Jegerlefmer, K. Schilling and R. Sommer (1994a): Nucl. 
Phys. B414, 815. 

C. Alexandrou, A. Borrelli, S. Giisken, F. Jegerlehner, K. Schilling, G. Siegert and R. 

Sommer (1994b): Phys. Lett. B337, 340. 
M. G. Alford, T. R. Klassen and G. P. Lepage (1997): Nucl. Phys. B (Proc. Suppl. 53), 

861. 



The Heavy Hadron Spectrum 



61 



A. Ali Khan in : Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B (Proc. 
Suppl.) to appear. 

A. Ali Khan, C. Davies, S. Collins, J. Sloan and J. Shigemitsu (1996): Phys. Rev. D53, 
6433. 

S. Aoki et al (1997): Nucl. Phys. B (Proc. Suppl. 53), 355. 

T. Appelquist, M. Dine and I. Muzinich (1978): Phys. Rev. D17, 2074. 

G. Bali, K. Schilling and A. Wachter (1997a): hep-lat/9703019, Phys. Rev. D (in press). 

G. Bali, K. Schilling and A. Wachter (1997b): Phys. Rev. D55, 5309. 

G. Bali and K. Schilling (1992): Phys. Rev. D46, 2636; ibid~D47, (1993) 661. 

A. Barchielli, N. Brambilla and G. Prosperi (1990): Nuovo Cimento 103A, no. 1, 59. 

G. Bodwin, E. Braaten and P. Lepage (1995): Phys. Rev. D51, 1125. 

P. A. Boyle (1997a): in Nucl. Phys. B (Proc. Suppl. 53), 398. 

P. A. Boyle (1997b): in Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B 

(Proc. Suppl.) to appear. 
W. Biichmuller and S.-H. H. Tye (1981): Phys. Rev. D24, 132. 

F. Butler, H. Chen, J. Sexton, A. Vaccarino and D. Weingarten (1994): Nucl. Phys. 
B430, 179. 

Y.-Q. Chen, Y. P. Kuang and R. J. Oakes (1995): Phys. Rev. D52, 264. 

F. Close (1979): An Introduction to Quarks and Partons, Academic Press. 

F. Close (1997): in Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B (Proc. 

Suppl.), to appear. 
S. Collins (1997): Nucl. Phys. B (Proc. Suppl. 53), 389. 

S. Collins (1997b): in Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B 

(Proc. Suppl.), to appear. 
S. Collins, R. G. Edwards, U. M. Heller and J. Sloan (1996a): Nucl. Phys. B (Proc. 

Suppl. 47), 455. 

S. Collins, U. Heller, J. Sloan, J. Shigemitsu, A. Ali Khan and C. Davies (1996b): Phys. 
Rev. D54, 5777. 

C. T. H. Davies (1997): in Proceedings of the International Workshop on Lattice QCD 
on Parallel Computers, Nucl. Phys. B (Proc. Suppl.), in press, hep-lat/9705039 and 
in Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B (Proc. Suppl.), to 
appear. 

C. T. H. Davies, K. Hornbostel, A. Langnau, G. P. Lepage, A. Lidsey, J. Shigemitsu 

and J. Sloan (1994a): Phys. Rev. D50, 6963. 
C. T. H. Davies, K. Hornbostel, A. Langnau, G. P. Lepage, A. Lidsey, C. Morningstar, 

J. Shigemitsu and J. Sloan (1994b): Phys. Rev. Lett. 73, 2654. 
C. T. H. Davies, K. Hornbostel, G. P. Lepage, A. Lidsey, J. Shigemitsu and J. Sloan 

(1995a): Phys. Lett. B345, 42. 
C. T. H. Davies, K. Hornbostel, G. P. Lepage, A. Lidsey, J. Shigemitsu and J. Sloan 

(1995b): Phys. Rev. D52, 6519. 
C. T. H. Davies, K. Hornbostel, G. P. Lepage, A. Lidsey, J. Shigemitsu and J. Sloan 

(1996) : Phys. Lett. B382, 131. 

C. T. H. Davies, K. Hornbostel, G. P. Lepage, P. McCallum, J. Shigemitsu and J. Sloan 

(1997) : Phys. Rev. D56, 2755. 

P. de Forcrand and J. Stack (1985): Phys. Rev. Lett. 55, 1254. 

DELPHI collaboration (M. Feindt) (1995): Invited talk at HADRON'95, CERN-PPE- 
95-139. 

DELPHI collaboration (1997a): Phys. Lett. B398, 207. 



62 Christine Davies 



DELPHI collaboration (1997b): Contribution to EPS-97 conference, pa-01,#452, avail- 
able at http://wwwcn.cern.ch/ pubxx/www/delsec/conferences/jerusalem. 

T. Draper, C. McNeile and C. Nenkov (1995): Nucl. Phys. B (Proc. Suppl. 42), 325. 

T. Draper and C. McNeile (1996): Nucl. Phys. B (Proc. Suppl. 47), 429. 

A. Duncan, E. Eichten, A. El-Khadra, J. Flynn, B. Hill and H. Thacker (1993): Nucl. 
Phys. (Proc. Suppl. 30), 433. 

A. Duncan, E. Eichten, J. Flynn, B. Hill, G. Hockney and H. Thacker (1995): Phys. 
Rev. D51, 5101. 

E852, D. R. Thompson et al (1997): Phys. Rev. Lett. 79, 1630. 

E. Eichten (1980): Phys. Rev. D22, 1819. 

E. Eichten, K. Gottfried, T. Kinoshita, J. Kogut, K. D. Lane and T.-M. Yan (1975): 

Phys. Rev. Lett. 34, 369. 
E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane and T. M. Yan (1978): Phys. Rev. 

D17, 3090 and ibid D21 (1980), 203. 
E. Eichten and F. Feinberg (1981): Phys. Rev. D23, 2724. 
E. Eichten and B. Hill (1990): Phys. Lett. B243, 427. 
E. Eichten, C. Hill and C. Quigg (1993): Phys. Rev. Lett. 71, 4116. 
E. Eichten and C. Quigg (1994): Phys. Rev. D49, 5845. 

A. X. El-Khadra and B. P. Mertens (1995): Nucl. Phys. B (Proc. Suppl. 42), 406. 
A. X. El-Khadra, A. S. Kronfeld and P. B. Mackenzie (1997): Phys. Rev. D55, 3933. 
A. Falk (1997): hep-ph/9707295. 

A. F. Falk, B. Grinstein and M. E. Luke (1991): Nucl. Phys. B 357, 185. 
S. S. Gershtein, V. V. Kiselev, A. K. Likhoded and A. V. Tkabladze (1995): Phys. Rev. 
D51, 3613. 

V. Gimenez, G. Martinelli and C. Sachrajda (1997): Nucl. Phys. B486, 227. 
M. Gorbahn, C. T. H. Davies, G. P. Lepage, J. Shigemitsu and J. Sloan (1997): in 
preparation. 

A. K. Grant, J. L. Rosner and E. Rynes (1993): Phys. Rev. D47, 1981. 

B. Grinstein and I. Rothstein (1996): Phys. Lett. B385, 265. 
D. Gromes (1977): Nucl. Phys. B131, 80. 

D. Gromes (1984): Z. Phys. C22, 265. 

D. Gromes (1988): Phys. Lett. B202, 262. 

H. Grosse and A. Martin (1980): Phys. Rep. 60C, 341. 

R. Gupta and T. Bhattacharya (1997): Phys. Rev. D55, 7203. 

G. Heatlie, C. T. Sachrajda, G. Martinelli, C. Pittori and G. C. Rossi (1991): Nucl. 
Phys. B352, 266. 

J. Hein (1997): in Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B (Proc. 
Suppl.) to appear. 

A. B. Henriques, B. H. Kellett and R. G. Moorhouse (1976): Phys. Lett. B64, 85. 

A. Huntley and C. Michael (1986): Nucl. Phys. B270, 123. 

A. Huntley and C. Michael (1987): Nucl. Phys. B286, 211. 

N. Isgur (1997): preprint JLAB-THY-97-26. 

N. Isgur and M. B. Wise (1991): Phys. Lett. B66, 1130. 

C. Itzykson and J.-B. Zuber (1980): Quantum Field Theory, McGraw Hill. 
A. S. Kronfeld (1997): Nucl. Phys. B (Proc. Suppl. 53), 401. 

K. J. Juge, J. Kuti and C. J. Morningstar (1997): in Proceedings of LAT97, Edinburgh, 

July 1997, Nucl. Phys. B (Proc. Suppl.), to appear. 
W. Kwong and J. L. Rosner (1988): Phys. Rev. D38, 279. 

R.Landua (1997): Review talk in Proceedings of ICHEP'96, Warsaw, World Scientific. 



The Heavy Hadron Spectrum 



63 



G. P. Lepage (1989): From Action to Answers, World Scientific. 
G. P. Lepage (1992): Nucl. Phys. B (Proc. Suppl. 26), 45. 

G. P. Lepage (1996): Lectures given at the 1996 Schladming Winter School, World 
Scientific. 

G. P. Lepage (1997): in Proceedings of Lattice QCD on Parallel Computers, Tsukuba, 

Nucl. Phys. B (Proc. Suppl.) in press, hep-lat/9707026. 
G. P. Lepage and P. B. Mackenzie (1993): Phys. Rev. D48, 2250. 

G. P. Lepage, L. Magnea, C. Nakhleh, U. Magnea and K. Hornbostel (1992): Phys. 
Rev. D46, 4052. 

M. Liischer, S. Sint, R. Sommer. P. Weisz and U. Wolff (1997): Nucl. Phys. B491, 323. 
T. Manke (1997): in Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B (Proc. 
Suppl.), to appear. 

T. Manke, I. T. Drummond, R. R. Horgan and H. P. Shanahan (1997): Phys. Lett. 

B408, 308. 
A. Martin (1980): Phys. Lett. B93, 338. 

A. Martin and J.-M. Richard (1987): Phys. Lett. B185, 426. 
C. Michael and P. E. L. Rakow (1985): Nucl. Phys. B256, 640. 

I. Montvay and G. Miinster (1994): Quantum Fields on a Lattice, Cambridge University 
Press. 

C. Morningstar (1994): Phys. Rev. D50, 5902. 
M. Neubert (1994): Phys. Rep. 245, 259 

S. Ono, A. Sanda and N. Tornqvist (1986): Phys. Rev. D34, 186. 
The Particle Data Group (1996), R. M. Barnett et at 

The Particle Data Group (1997): updated particle tables at http://pdg.lbl.gov/. Phys. 
Rev. D 54, 1. 

S. J. Pcrantonis and C. Michael (1990): Nucl. Phys. B347, 854. 

M. Peskin (1983): in 11th SLAC Summer Institute, SLAC Report PUB 3273. 

J. Peisa and C. Michael (1997): in Proceedings of LAT97, Edinburgh, July 1997, Nucl. 

Phys. B (Proc. Suppl.) to appear. 
C. Quigg (1997a): Physics Today, Vol. 50 no. 5, 20. 
C. Quigg (1997b): hep-ph/9707493. 

C. Quigg and J. L. Rosner (1977): Phys. Lett. 71B, 153. 
C. Quigg and J. L. Rosner (1979): Phys. Rep. 56C, 167. 
P. A. Rapidis et al (1977): Phys. Rev. Lett. 39, 526. 
J. Richardson (1979): Phys. Lett. B82, 272. 
J. Rosner (1995): hep-ph/9501291. 

H. J. Schnitzer (1975): Phys. Rev. Lett. 35, 1540. 

SESAM collaboration, U. Glassner, S. Giisken, H. Hoeber, T. Lippert, G. Ritzenhofer, 
K. Schilling, G. Siegert, A. Spitz and A. Wachter (1996): Phys. Lett. B383, 98. 

SESAM collaboration, N. Eicker, Th. Lippert, K. Schilling, A. Spitz, J. Fingberg, S. 
Giisken, H. Hoeber and J. Viehoff (1997): hep-lat/9709002. 

H. Shanahan (1997): in Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B 
(Proc. Suppl.), to appear. 

B. Sheikholeslami and R. Wohlert (1985): Nucl. Phys. B259, 572. 

J. Simone et al (1997): in Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B 

(Proc. Suppl.), to appear. 
R. Sommer (1994): Nucl. Phys. B411 839. 
R. Sommer (1997): these Proceedings. 



64 Christine Davies 



A. Spitz (1997) in: Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B (Proc. 
Suppl.), to appear. 

B. A. Thacker and G. P. Lepage (1991): Phys. Rev. D43, 196. 

H. B. Thacker, E. Eichten and J. C. Sexton (1988): Nucl. Phys. B (Proc. Suppl. 4), 
234. 

H. Trottier (1997a): Phys. Rev. D55, 6844. 

H. Trottier (1997b) in: Proceedings of LAT97, Edinburgh, July 1997, Nucl. Phys. B 

(Proc. Suppl.), to appear. 
UKQCD collaboration, S. P. Booth et al (1992a): Phys. Lett. B284, 377. 
UKQCD collaboration, S. P. Booth, D. S. Henty, A. Hulsebos, A. C. Irving, C. Michael 

and P. W. Stephenson (1992b): Phys. Lett. B294, 385. 
UKQCD collaboration, A. K. Ewing et al (1996a): Phys Rev. D 54, 3526. 
UKQCD collaboration, K. C. Bowler et al (1996b): Phys Rev. D 54, 3619. 
UKQCD collaboration, presented by R. D. Kenway (1997): Nucl. Phys. B (Proc. Suppl. 

53), 209. 

WA89 (R. Werding) (1995) : Proceedings of ICHEP'94, Glasgow, IOP Publishing. 
D. Weingarten (1997): these Proceedings. 



